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ABSTRACT 

Middle  ultraviolet  spectra  from  1800  to  3400A  are 
analyzed.  These  spectra  were  obtained  from  the  March  1990 
rocket  flight  of  the  NPS  MUSTANG  instrument  over  the  altitudes 
105km  to  315km.  The  data  were  compared  with  computer 
generated  synthetic  spectra.  A  least  squares  fitting 
procedure  was  developed  for  this  purpose.  Each  data  point  was 
weighted  using  the  standard  deviation  of  the  means.  Synthetic 
spectra  were  generated  for  the  following  emissions:  N2 
Vegard-Kaplan;  N2  Lyman-Birge-Hopf ield;  NO  y,  8,  and  e;  01 
2912k,  Oil  2470A;  and  Nil  2143A.  Altitude  profiles  for  the 
emissions  were  obtained.  Tentative  identification  was  made  of 
the  OIII  2853A  emission.  A  comparison  of  VK  and  LBH  profiles 
demonstrates  the  process  of  N2  A-state  quenching  by  atomic 
oxygen. 
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I.   INTRODUCTION 

The  composition  of  the  earth's  atmosphere  is  important 
to  many  areas  of  science,  such  as  weather  prediction, 
satellite  operations,  and  communications.  Properties  of  the 
upper  atmospheric  region  known  as  the  ionosphere  are  of 
specific  interest  to  the  military,  as  they  strongly  influence 
such  systems  as  high  frequency  (HF)  radio  communications, 
over-the-horizon  radar  (OTH) ,  ballistic  missile  early  warning 
(SDI) ,  and  the  Ground  Wave  Emergency  Network  (GWEN) .  The 
electron  density  profile,  a  plot  of  electron  density  versus 
altitude  for  a  given  time  and  location,  is  needed  for 
effective  use  of  the  ionosphere.  Electron  density  is  a  factor 
that  varies  not  only  with  altitude,  but  with  time  and 
geographical  location  as  well.  In  1986,  the  Joint  Chiefs  of 
Staff  prioritized  the  determination  of  ionospheric  electron 
density  as  fifth  on  a  list  of  43  critical  geophysical 
parameters  for  inclusion  in  a  defense  satellite  (The  Joint 
Chiefs  of  Staff, 1986). 

Currently,  electron  density  profiles  are  determined  using 
ground-based  ionosonde  stations.  This  is  very  effective 
locally,  but  is  inherently  inaccurate  for  global  forecasting 
since  there  are  only  about  20  stations  worldwide.  A 
satellite-based  method  of  forecasting  the  ionosphere  would  be 


globally  effective,  but  to  place  an  ionosonde  in  orbit  is 
unfeasible  due  to  size  and  power  limitations.  One  alternative 
to  ionosondes  is  to  infer  electron  density  profiles  from 
spectrographic  measurements  of  the  atmospheric  airglow.  An 
operational  spectrograph  is  small  and  relatively  low-powered; 
it  can  easily  be  included  in  a  defense  satellite  package.  The 
Naval  Research  Lab  and  the  Naval  Postgraduate  School  have  a 
joint  spectrographic  project  scheduled  to  be  included  in  a 
satellite  to  be  launched  in  1995.  The  NRL  project  is  a  high 
resolution  airglow  and  aurora  spectrograph  (HIRAAS)  and  the 
accompanying  NPS  effort  is  the  middle  ultraviolet  spectrograph 
(MUSTANG) .  The  inclusion  of  the  duo  represents  a  significant 
improvement  in  the  global  determination  of  electron  density 
profiles  over  ground-based  ionosondes. 

The  inference  of  electron  densities  from  emission 
observations  by  the  HIRAAS /MUSTANG  is  performed  by  relating 
measurements  of  natural  ultraviolet  emissions  to  the  neutral 
species'  densities.  A  rocket-launched  HIRAAS/MUSTANG 
experiment  was  first  tested  in  March  1990.  The  data  from  the 
launch  is  the  topic  of  this  thesis. 

A.   THESIS  OBJECTIVES 

This  thesis  focuses  on  the  analysis  of  data  obtained  from 
the  March  1990  flight.  The  analysis  extends  methods  used 
previously  on  subsets  of  the  data  to  the  entire  observed 
range,   180oA  to  340oA.    All  known  synthetic  models  with 


emissions  in  this  range  are  used  in  the  fitting.  The  result 
of  the  fitting  algorithm  is  a  set  of  relative  scale  factors 
for  the  various  spectra  by  altitude.  The  relative  scale 
factors  are  converted  into  absolute  terms  and  plotted  against 
altitude.  These  absolute  plots  of  intensity  and  column 
density  can  be  used  to  infer  an  electron  density  profile  (see 
e.g.,  Meier, 1991;  McCoy, 1985).  A  secondary  goal  is  to 
determine  the  feasibility  of  an  approximation  for  the  N2 
Second  Positive  spectrum  given  the  unavailability  of  a 
synthetic  spectrum.  The  determination  of  the  Second  Positive 
spectrum  would  allow  an  estimation  of  photoelectron  flux, 
since  it  is  produced  by  transitions  from  a  state  which  can 
only  be  populated  by  collisions. 

B.   THESIS  OUTLINE 

The  thesis  is  divided  into  six  chapters.  Chapter  II  gives 
general  background  information  on  the  descriptions  and 
designations  of  the  atmosphere  and  existing  emissions.  A 
basic  explanation  of  atomic  and  molecular  spectra  is  also 
summarized  in  this  chapter.  Chapter  III  addresses  the 
specific  emission  features  in  the  wavelength  range  of 
interest.  Chapter  IV  presents  a  summary  of  the  experiment  and 
the  subsequent  data  reorganization.  Chapter  V  gives  a 
complete  treatment  of  analysis  of  the  data  and  results  of  the 
study.  Chapter  VI  concludes  the  paper  with  a  summary  of 
findings  and  suggestions  for  future  study.   The  Appendices 


contain  the  computer  programs  which  generated  the  synthetic 
spectra  and  the  fitting  programs,  as  well  as  plots  of  the 
final  fits  to  the  data. 


II.   BACKGROUND 

A.   THE  ATMOSPHERE 

1.   General  Description 

The  Earth's  atmosphere  is  designated  by  horizontal  layers, 
and  these  divisions  are  determined  by  one  of  two  primary 
methods.  One  method  designates  the  layers  based  upon 
temperature  characteristics;  the  other  uses  density  to 
categorize  the  strata.  The  layers  and  a  typical  temperature 
profile  are  shown  in  Figure  2-1. 
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Figure  2-1  Temperature  profile  and  corresponding  atmospheric 
layers. 


With  respect  to  the  temperature  designations,  the 
troposphere  is  the  region  closest  to  the  Earth's  surface  where 
temperatures  decrease  with  altitude,  typically  from  about  3  00K 
at  the  surface  to  220K  at  10  km.  The  point  at  which 
temperature  begins  to  rise  again  is  called  the  tropopause  and 
it  indicates  the  start  of  the  next  region,  the  stratosphere. 
The  stratosphere  is  a  relatively  warm  layer  due  to  absorption 
of  solar  radiation (primarily  2000-3000A)  by  ozone.  The  top  of 
the  stratosphere  is  marked  by  the  stratopause.  At  this  point 
the  mesosphere  begins.  The  temperature  again  begins  to 
decrease  with  altitude  until  a  low  of  about  100K  occurs  at 
around  85  kilometers.  This  is  the  mesopause ,  which  marks  the 
beginning  of  the  next  layer,  the  thermos  phere .  In  the 
thermosphere,  temperature  rises  steadily  with  increasing 
altitude  until  reaching  an  eguilibrium  temperature  between 
1000K  and  2000K,  depending  on  solar  activity. 

Studying  the  atmosphere  in  terms  of  density  gives  rise  to 
a  different  set  of  atmospheric  layers.  Under  this  convention, 
the  lowest  80  kilometer  layer  of  the  atmosphere  is  called  the 
homosphere .  In  this  region,  vertical  mixing  of  the  air 
currents  provides  a  uniform  relative  density  ratio  between 
various  species.  The  primary  components  are  nitrogen,  oxygen, 
and  argon  in  a  ratio  of  78:21:0.9.  The  remaining  0.1  percent 
is  primarily  helium,  hydrogen,  and  carbon  dioxide.  The  20  km 
from  80  to  100  is  designated  the  harosphere ,  or  tropopause, 
and  it  is  a  region  of  transition  above  which  very  little 


vertical  mixing  occurs.  Above  the  barosphere,  from  100  km  to 
around  1000  km,  lies  the  heterosphere .  Due  to  the  absence  of 
mixing  in  this  region,  the  ratio  of  constituents  varies  with 
altitude.  At  extremely  high  altitudes  of  1000  kilometers  or 
so,  the  atmosphere  becomes  so  tenuous  that  particle 
trajectories  are  essentially  parabolic  orbits.  This  "outer 
atmosphere"  is  called  the  exosphere. 

The  data  analysis  in  this  paper  is  concentrated  in  the 
lower  heterosphere,  the  first  2  00  km  of  the  thermosphere, 
where  there  is  an  obvious  positive  temperature  gradient  and 
the  convective  mixing  is  minimal.  This  region  is  also  termed 
the  ionosphere. 

The  ionosphere  is  so  termed  due  to  the  ionization  of  the 
atmosphere  by  absorbed  solar  radiation.  The  actual  height  and 
thickness  of  the  four  layers  (D,  E,  Fl,  F2)  vary  with  solar 
activity  and  atmospheric  composition.  The  D  region  is 
typically  from  70  km  to  90  km  during  the  day;  the  E  region  is 
typically  from  90  to  130  km;  Fl  generally  exists  from  130  to 
160  km  during  the  day;  and  F2  is  above  160  km.  Each  of  these 
regions  is  formed  by  ionization  due  to  solar  radiation  at 
different  wavelengths.  For  example,  longer  wavelength 
radiation  can  penetrate  deeper  into  the  atmosphere.  Since 
solar  activity  changes  daily,  seasonally,  and  yearly,  the 
intensity  of  sunlight  at  a  particular  wavelength  and  hence  the 
depth  of  penetration  is  variable.  Ionization  depends  on  the 
composition  of  the  neutral  atmosphere  as  well.    Combining 


these  two  factors  of  composition  and  solar  activity,  one  can 
see  that  the  flux  of  radiation  is  greatest  at  high  altitude 
but  density  is  lowest  and  vice  versa.  The  maximum  ionization 
occurs  somewhere  between  these  two  extremes,  and  the  result  is 
an  overall  layer  known  as  a  Chapman  layer.  The  composition 
and  formation  of  the  ionosphere  is  quite  complex  and  is  not 
unique  to  Earth.  An  excellent  discussion  of  planetary 
ionospheres  is  written  by  Chamberlain  (1978) . 

2 .   Atmospheric  Airglow 

Airglow  is  the  naturally  occurring  radiation  in  the 
Earth's  atmosphere.  The  daytime  emissions,  called  dayglow, 
are  primarily  driven  by  photoexcitation.  Also  contributing  to 
dayglow  are  emissions  which  result  from  photoelectron  impacts 
and,  to  a  lesser  degree,  photochemical  reactions.  During 
daylight  hours  photochemical  reactions  are  negligible  compared 
to  the  abundance  of  emissions  due  to  sunlight,  but  the  same 
reactions  are  significant  at  night  since  there  are  obviously 
few  direct  solar  photons  to  stimulate  transitions.  Because  the 
data  analyzed  in  this  work  were  obtained  during  the  day, 
photochemical  reactions  are  not  specifically  addressed  in  this 
work. 

The  primary  source  of  dayglow  emissions  is 
photoexcitation.  Photoexcitation  is  the  absorption  of  a 
photon  as  a  result  of  inelastic  collision  between  a  photon  and 
an  atom  or  molecule  to  raise  the  particle  to  an  excited  state. 


When  the  photon  collides  with  the  particle,  the  photon  is 
absorbed  and  the  particle  jumps  to  a  higher  energy  level. 
When  the  particle  drops  back  to  a  lower  energy  state,  it  emits 
a  photon  with  energy  equal  to  the  transition.  This  process 
will  be  further  described  in  the  next  section  on  atomic  and 
molecular  emissions. 

Although  not  as  important  to  dayglow  as  photoexcitation, 
photoelectron  collisions  are  significant  because  they  can 
produce  transitions  which  are  otherwise  forbidden. 
Observations  of  forbidden  transitions  can  be  used  to  find  the 
photoelectron  flux,  since  only  photoelectron  collisions  could 
excite  the  particle  to  the  upper  state.  Given  the  importance 
of  whether  a  transition  is  allowed  or  forbidden  to 
calculations  of  photoelectron  flux,  a  knowledge  of  basic 
selection  rules  is  desirable  for  the  analysis  of  dayglow 
emissions. 

It  is  also  worthy  to  note  that  some  of  the  forbidden 
transitions  which  are  prominent  in  the  low  density  unconfined 
gas  of  the  thermosphere  cannot  be  produced  in  earthbound  lab 
experiments  (Rees,1989).  Such  emissions,  seen  in  dayglow  data 
collected  by  space-borne  equipment,  do  not  match  any  of  the 
spectra  which  have  been  observed  using  lab  or  ground-collected 
data.  These  emissions  are  identified  with  transitions  of 
known  elements,  by  wavelength  matching.  One  needs  to  be 
familiar  with  spectral  nomenclature  and  selection  rules  in 
order  to  understand  the  explanations  of  these  identifications. 


B.   ATOMIC  AND  DIATOMIC  SPECTRA 
1.   Atomic  Emissions 

The  primary  atom-photon  interactions  are  shown  in  Figure 
2-2.  In  elastic  scattering,  no  energy  transfer  takes  place 
between  the  photon  and  the  atom  during  collision,  although  the 
photon  may  change  directions.  In  inelastic  scattering,  the 
incident  photon  may  not  only  change  direction,  but  may  impart 
energy  to  the  atom.  In  this  case,  the  atom  is  excited  to  a 
level  AE  above  its  original  state  and  the  deflected  photon 
will  have  an  energy  AE  less  than  its  incident  energy.  A 
photon  with  energy  AE  will  be  emitted  when  the  atom  relaxes  to 
the  initial  state.  Resonance  scattering  occurs  when  the 
incident  photon  is  absorbed,  raising  the  atom  to  an  excited 


a)  Elastic  scattering 


b)  Inelastic  scattering 


c)  Resonance  radiation 


d)  Fluorescence 


e)  Collisional  De-excitation 


Figure  2-2   Atomic  processes  induced  by  incident  radiation 

(Brehm  &  Mullin,  1989) . 
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state.  A  photon  with  wavelength  equal  to  the  incident  photon 
is  emitted  when  the  atom  falls  back  to  its  initial  state. 
Fluorescence  is  the  phenomenon  which  occurs  when  the  excited 
atom  undergoes  several  downward  transitions  in  returning  to 
the  initial  state,  giving  a  cascade  of  photons.  An  atom  can 
also  return  to  the  ground  state  by  colliding  with  another 
particle,  thus  returning  to  the  ground  state  without  emitting 
a  photon.  This  is  called  collisional  de-excitation  or 
quenching. 

Similar  to  photon-atom  collisions,  electron-atom 
collisions  may  be  elastic  or  inelastic.  In  an  electron-atom 
collision,  the  scattered  electron  loses  kinetic  energy,  rather 
than  actually  being  absorbed,  when  the  atom  is  excited  to  a 
higher  energy  level.  Photon  emissions  occur  when  the  atom 
relaxes  to  the  ground  state. 

2.   Molecular  Emissions 

Molecular  spectra  may  occur  as  the  result  of  electronic 
state  changes,  but  they  may  also  result  from  purely 
vibrational  and  rotational  state  changes.  Each  of  these  types 
of  transitions  has  typical  energy  gaps,  with  electronic  being 
on  the  order  of  10  eV,  vibrational  of  .1  eV,  and  rotational  of 
.001  eV.  These  three  types  of  transitions  are  illustrated  in 
Figure  2-3.  As  seen  in  the  figure,  for  a  given  electronic 
energy  level  there  are  many  vibrational  energy  levels. 
Similarly,  for  a  given  change  in  vibrational  energy  there  are 
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Figure  2-3   Molecular  absorption  transitions: 
Transition  I  is  pure  rotational 
Transition  II,  vibrational  change  A^=l 
Transition  III,  electronic  change  of  state 
(Brehm  &  Mullin,  1989) . 

many  possible  initial  and  final  rotational  levels.  As  a 
result  of  the  many  vibrational  and  rotational  levels  which  can 
be  populated  for  any  given  electronic  transition,  the  emission 
feature  is  composed  of  many  closely  spaced  lines.  Figure  2-4 
shows  a  molecular  spectrum,  exhibiting  a  distribution  of 
emissions.  When  the  separation  of  component  lines  is  smaller 
than  instrument  resolution,  the  transition  may  appear  as  a 
hump,  as  shown  in  Figure  2.5.  This  is  the  case  for  the  data 
studied  in  this  paper. 
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The   study   of   quantum 

mechanics     explores 

transitional    possibilities, 

and    develops    a    set    of 

guidelines  for  those  emissions 

one  expects  to  observe   These 

so-called    selection    rules 

require    that    rotational 

Figure    2-4        Graphical 

transitions   only   occur   for  representation   of   molecular 

spectrum. 

aj=±1,  where  j  is  the  total 

quantum    number 

( j=0, 1, 2, 3 , . . . ) .  Vibrational 

transitions  occur  within  the 

confines  of  selection  rules 

aj=±1  and  av=±1,  where  v  is  Figure  2-5  Instrument  response 

to  molecular  transition. 

the    vibrational    quantum 

number  (v=0,l, . . . j) .  These  rules  are  not  strictly  adhered  to, 
as  they  are  based  on  electronic  dipole  behavior,  and  some 
molecules  exhibit  quadrupole  and  magnetic  dipole 
characteristics.  There  are  also  cases  where,  though  only 
slightly  probable,  emissions  occur  when  av>±1.  Electronic 
transition  selection  rules  are  complex,  and  a  full  treatment 
is  given  to  the  development  and  explanation  of  them  by  Green 
and  Wyatt  (1965).   A  summary  of  these  is  shown  in  Table  2-1. 
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TABLE   2-1   SUMMARY   OF   SELECTION   RULES   FOR   MOLECULAR 
TRANSITIONS  (Green  &  Wyatt,  1965) . 


Electric  Dipole  (El) 

Magnetic  Dipole  (Ml) 

Electric  Quadrupole  (E2) 

Total  inversion 
Electronic 
inversion 

+  ♦*  +    —  *-*  — 
g<*  g       u  <->  u 

g  <-+  g       u*-hl 

AJ 
A5(a),(b) 

0,±1 
0 

0*,±l(0«+»0) 
0 

Q,  ±1,  ±2t 

0 

A2(a} 

0 

SeeAA 

0 

AA(a) 

0,  ±1 

ilifAS  =  0, 
Oif  AS  =  ±1 

0,  ±1,  ±2 

AQ(a) 

Qfc±l 

±1 

0,  ±1,  ±2 

AA  (b) 

0*,±[ 

0,±1 

0,  ±1,  ±2 

Atf(b) 

oui 

o,±i§ 

0,  ±1§,  ±2 

•NotO  -0.      f  Not  0-  0,0-1,  i- 1 

J  If  0  -  0  then  1/  -  0.       §  Not  for  2  -  I  transitions. 


Derivation  of  these  rules  is  not  essential  to  this  data 
analysis,  but  an  understanding  of  the  corresponding  molecular 
nomenclature  is  necessary  in  order  to  understand  the 
discussion. 

The  component  of  total  angular  momentum  which  lies  along 
the  internuclear  axis  of  a  molecule  is  designated  A,  and  is 
used  to  label  molecular  electronic  state.  Each  state  is  given 
a  capital  Greek  letter  depending  on  its  value  of  A.  These 
designations  are  shown  in  Table  2-2.  The  spin  of  the  molecule 
is  designated  S  (which  should  not  be  confused  with  the  E 
symbol  for  the  A=0  state) .  It  can  take  on  values  from  -S  to  S, 
where  S  is  the  resultant  of  individual  electron  spins  in  the 
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TABLE     2-2     MOLECULAR 
ELECTRONIC  STATES. 


molecule.  Multiplicity ,  which  is 
2S+1,  is  written  as  a  left 
superscript  on  the  state  letter. 
The  ground  state  of  a  diatomic 
molecule  is  labelled  with  a 
leading  X.  Successive  excited 
states  of  the  same  multiplicity 
are  labelled  A,B,C,  and  so  forth. 
Excited  states  of  different 
multiplicity   are   alphabetically 


DESIGNATIONS  OF 
MOLECULAR  STATES 

Greek 
Letter 

A  of 
Molecule 

£ 

0 

n 

±1 

A 

±2 

$ 

±3 

designated  by  small  letters.  If  these  labels  are  strictly 
followed,  then  in  general  the  selection  rules  allow 
transitions  between  like-labelled  states  (e.g. 
capital^capital)  .  As  with  any  rule,  there  are  exceptions. 
For  instance,  N2  diagrams  have  capital  and  small  letters 
assigned  conversely. 
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III.   DAYGLOW  EMISSION  FEATURES 


A.   MOLECULAR  NITROGEN 

Figure  3-1  is  an  N2  energy  level  diagram  including  the 
three  transition  bands  of  interest  in  the  1800-3400A  range. 
The  states  are  indicated  by 
molecular  notation,   and  the 
transitions  are  identified  by 
name. 

1.   Vegard-Kaplan 

The  Vegard-Kaplan  system 
is  so  designated  because 
Vegard  discovered  the  system 
in  1932,  and  Kaplan  determined 
the  specific  wavelengths 
involved  a  few  years  later. 
The  emissions  range  from  1250A 
to  5325A  and  represent  the 
A32,.+->X1Z,.+  transition.  In  this 


Figure  3-1 


energy  level 


transition,     the    total  ^^am  (McEwan,  1975). 
transformation  is  +-*•+,  the  electronic  inversion  is  u->g,  AE=-1, 
and  AA=0.   From  the  molecular  selection  rules  summarized  in 
Table  2-1,  this  is  a  forbidden  transition,  yet  it  has  been 
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observed  as  an  electric  dipole  transition.    A  thorough 
discussion  of  the  band  system  is  given  in  Meier  (1990). 

A-state  nitrogen  molecules  are  produced  by  cascading 
downward  from  the  B  and  C  states  or  as  a  result  of 
photoelectron  collisions.  Because  the  relationship  of 
emission  bands  to  photoelectron  flux  is  complicated  by 
cascading  N2  molecules  and  quenching  effects,  the  V-K  band  is 
not  generally  used  to  infer  photoelectron  flux. 

2.  Lyman-Birge-Hopf ield 

The  Lyman-Birge-Hopf ield  system  is  the  name  given  to  the 
a1IIg->X1Z  +  transition  of  N2.  The  emission  bands  arising  from 
this  transition  extend  from  1160A  to  302oA  and  are  due  to  both 
electric  quadrupole  and  magnetic  dipole  transitions.  Unlike 
the  V-K  system,  the  upper  state  of  the  LBH  system  can  only  be 
populated  by  photoelectron  collision  excitation.  This  system 
can  therefore  be  used  to  directly  infer  the  photoelectron 
flux.  The  LBH  bands  found  in  the  data  were  particularly 
significant  in  the  wavelength  region  from  I8O0A  to  2400A.  An 
intensive  investigation  of  the  LBH  system  as  it  pertains  to 
the  MUSTANG  data  was  performed  by  Mack  (1991) . 

3.  Second  Positive 

The  Nitrogen  Second  Positive  system  (C3IIu->B2ng)  extends 
from  2680A  to  5460A.  It  is  a  particularly  important 
transition  because  the  upper  level  is  not  directly  populated 
through  photoexcitation,  but  only  by  collision.    Thus  a 
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measurement  of  the  Second  Positive  intensity  profile  can  be 
used  to  infer  the  photoelectron  flux.  This  result  can  be 
compared  with  those  obtained  from  the  LBH  bands. 
Unfortunately,  there  is  no  synthetic  spectrum  model  developed 
for  the  Second  Positive.  It  was  originally  intended  to 
determine  the  N2  Second  Positive  spectrum  by  subtracting  known 
contributions  from  the  data  and  assuming  the  remaining  data 
was  Second  Positive.  This  means  of  estimating  a  spectrum  was 
found  to  be  impractical  due  to  inaccuracies  in  the  data  caused 
by  data  drops. 

B.   ATOMIC  NITROGEN 

Figure  3-2  shows  the  transitions  which  have  been  observed 
in  the  airglow  due  to  N+.  Two  transitions  are  in  the 
wavelength  range  of  this  study.  The  3070A  emission  falls  in 
the  same  range  of  wavelengths  as  the  N2  Second  Positive 
spectrum,  so  it  could  not  be  definitively  seen  in  the  data. 
The  emission  shown  at  2143A  was  considered  in  the  synthetic 
fit  to  the  MUSTANG  data.  It  is  actually  a  doublet  at  2143.  55A 
and  2139. 68A  with  a  ratio  of  100:58,  respectively.  This 
intensity  ratio  was  determined  by  Bucsela  and  Sharp  (1989) 
from  dayglow  observation.  Theoretical  and  laboratory 
estimates  of  relative  line  strength  range  from  100:60  to 
100:23.  Since  the  data  source  for  the  Bucsela  and  Sharp  ratio 
was  similar  to  the  experimental  data  being  analyzed,  the 
100:58  relation  was  used  in  the  synthetic  model. 
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Figure  3-2   Energy  level  diagram  for  N+  (Rees,  1989)  . 

C.   ATOMIC  OXYGEN 

Figures  3-3  and  3-4  diagram  the  atomic  oxygen  emissions 
which  have  been  viewed  in  the  airglow.  Two  of  these  emissions 
fall  in  the  wavelength  range  explored  by  this  analysis. 

The  neutral  oxygen  atom  emits  a  photon  with  wavelength 
2972A  when  it  falls  from  the  1S  state  to  the  2P  ground  state. 
This  transition  is  depicted  in  Figure  3-3.  Singly  ionized 
oxygen,  shown  in  Figure  3-4,  produces  a  line  at  2470A  when  it 
falls  from  the  2P  state  to  the  4S  ground  state. 
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Figure  3-3   Energy  level  diagram  for  neutral  atomic  oxygen 
(Rees,  1989) . 

D.   NITRIC  OXIDE 

Nitric  oxide  emissions,  shown  in  Figure  3-5,  dominate  much 
of  the  dayglow  data.  There  are  four  primary  bands  in  the 
range  180oA  -  3400A,  7,  S ,  (3 ,  and  e.  The  &  bands  are  several 
orders  of  magnitude  less  intense  than  the  other  emissions,  so 
they  have  been  neglected. 

1.   Delta 

The  C2II->X2n  transition  band  system,  designated  S,  ranges 
from  162lA  to  4288A  (Barth, 1965) .  Using  standard  notation, 
the  upper  electronic  state's  vibrational  level  is  designated 
v'  ,  and  the  lower  state  vibrational  level,  v" .  The  «5 
transition  band  includes  vibrational  levels  corresponding  to 
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v'=0, 1,2,3,4  and  v"=0 , 1 , . . . 23 .  Nitric  oxide  will  dissociate 
at  v'>0  for  C2I1,  so  v'  is  limited  to  v'=0  and  v"<14  in 
creating  the  synthetic  S    spectrum. 

2 .   Gamma 

The  7  bands  of  nitric  oxide  extend  from  1873A  to  6126A. 
The  emissions  arise  from  the  A2E+-+X2II  transition,  and  the 
ground  state  has  v"=0,...23  as  previously  stated.  The 
dissociation  energy  for  nitric  oxide  is  lower  than  the  energy 
of  the  v'=4  level  in  the  A  state,  so  only  transitions  arising 
from  v'=0, 1,2,3  contribute  to  the  synthetic  y   spectrum. 
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Figure  3-4   Energy  level  diagram  for  0+  (Rees,  1989) . 
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3.   Epsilon 

The  nitric  oxide  e  system 
is  the  result  of  molecules 
transitioning  from  D2E+-+X2II. 
The  wavelengths  of  emissions 
range  from  I8O0A  to  2385A.  D- 
state  molecules  populate 
vibrational  levels  v'=0,l,2. 


Figure  3-5    Transitions  for 
nitric  oxide  (McEwan,  1975) . 
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IV.   THE  EXPERIMENT 

A.  Introduction 

The  middle  ultraviolet  spectra  analyzed  in  this  thesis 
were  obtained  from  the  NPS  MUSTANG  spectrograph  flown  on  a 
NASA  sounding  rocket  launched  March  30,  1990,  from  White  Sands 
Missile  Range,  New  Mexico.  The  spectrograph  was  launched  on 
a  Terrier  boosted  Black  Brant,  which  carried  the  pay load  to  an 
apogee  altitude  of  32  0  km.  The  spectrograph  was  looking  in 
the  anti-solar  direction  with  an  observation  zenith  angle  of 
90°.  From  320  to  180  km  on  the  downleg  no  useful  data  were 
obtained  due  to  the  pay load  attitude. 

B.  Instrument  Description 

The  MUSTANG  (Middle  Ultraviolet  Spectrograph)  is  a 
modified  1/8  m  Ebert-Fastie  spectrograph  with  a  1/8  m 
telescope.  The  entrance  slit  is  5mm  x  140  /xm  and  a  loA 
wavelength  resolution  is  provided  by  a  1200  g/mm  grating 
ruling.  The  detector  assembly  consists  of  an  ultraviolet-to- 
visible  image  intensifier  using  a  quartz  input  window,  a  CsTe 
photocathode, and  a  fiber  optic  output  window.  The  detector  is 
an  array  of  512  photodiodes  giving  voltage  outputs 
proportional  to  intensity  using  a  50  msec  integration  period. 
This  voltage  output  is  fed  serially  to  an  A/D  converter  and 
then   to   rocket   telemetry   via   a   FIFO   shift   register. 
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Calibration/sensitivity  curves  were  obtained  using  a 
deuterium  lamp  for  the  wavelength  range  1800-2800  A  and  a  FEL 
tungsten  lamp  for  2400-3400A.  A  more  precise  calibration  of 
the  1800-2000A  range,  which  is  affected  by  02  absorption,  was 
performed  in  a  vacuum  chamber  (Mack, 1991). 

C.   Data  Collection 

Over  8000  spectra  were  collected  by  MUSTANG.  These 
spectra  were  then  averaged  into  2  5  altitude  bins,  disregarding 
spectra  which  were  contaminated  by  the  instrument's 
electronics.  Each  of  the  25  resultant  spectra  covers  a  10  km 
layer  and  is  identified  by  its  central  altitude. 

Because  of  a  timing  difference  between  clocking  and 
reading  functions,  only  480  of  the  512  photodiodes'  responses 
were  recorded.  It  was  first  determined  that  every  seventeenth 
word  was  dropped,  and  these  blanks  were  inserted  as  averages 
of  the  two  adjoining  data  values  (Clayton, 1990) .  This  still 
did  not  produce  a  suitable  wavelength  match  with  known 
emissions,  so  a  mathematical  spline  was  performed  on  the  data 
in  order  to  match  a  uniform  wavelength  grid  (Mack, 1991). 
Although  this  provided  an  excellent  fit  to  assumed  models,  it 
had  no  physical  justification.  Laboratory  testing  was 
performed  by  Quint  (1991)  to  determine  which  values  were  lost 
in  the  timing  error.  His  results  accounted  for  the  32 
expected  data  dropouts,  but  two  additional  dropouts  remained 
unexplained.    These  two  shifts  were  inserted  at  positions 
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which  were  consistent  with  the  findings  of  Quint  and  also 
provided  the  best  agreement  with  three  previous  studies.  It 
should  be  pointed  out  that  this  was  not  confirmed  with  the 
electronics.  The  corrected  spectra  are  delineated  in  Appendix 
A,  noting  each  dropped  pixel.  At  these  spots  an  average  of 
the  adjacent  values  was  inserted  into  the  original  480  data 
points,  to  obtain  an  array  of  512  values. 

The  adjusted  spectra  were  then  divided  by  instrument 
sensitivity.    The  first  six 


values  were  set  to  zero,  since 

the  uncertainty  in  sensitivity 

of   the   instrument   at   the 

endpoint  is  extremely  large. 

It  was  found  that  one  further 

correction  had  to  be  made  to 

the  data.   In  both  pre-launch 

and  post-launch  high  voltage 

dark    count    spectra,    two 

voltage  spikes  existed.  These 

hot  spots  occurred  at  points 

corresponding   to   2850A  and 

Figure   4-1     MUSTANG   "hot 

3373A,  as  shown  in  Figure  4-1  pixels';  plotted  as  intensity 

(in  R/A)  versus  wavelength. 

These  dark  count  hot  spots 

were  converted  to  R/A  and  then  subtracted  from  the  averaged 

spectra  before  analysis. 
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V.   DATA  ANALYSIS  &  RESULTS 

A.   TECHNIQUE 

1.   Introduction 

a.    Fitting  Method 

Synthetic  spectra  were  fit  to  the  data  using  a  least 
squares  fit  process  (hereafter  referred  to  as  GRIDLS) .  It  was 
developed  by  Bevington  (1969)  and  applied  in  limited  portions 
of  the  wavelength  range  by  Clayton  (1990) ,  Andersen  (1990) , 
and  Mack  (1991)  .  This  method  is  an  extension  of  a  basic 
least-squares  regression,  minimizing  error  between  the  fit  and 
the  data  as  measured  by  the  variable  x2  with  weighting  for 
confidence  in  the  data.  The  program  fits  a  linear  combination 
of  synthetic  spectra  to  the  data,  with  each  synthetic  spectrum 
having  an  independent  coefficient,  or  scale  factor.  The 
program  is  adaptable  for  any  number  of  desired  parameters. 
The  final  analysis  used  nine  parameters,  which  will  be 
discussed  in  the  Results  section.  The  fitting  program 
(GRIDLS)  and  its  accompanying  subprograms  are  in  Appendix  B. 

The  output  of  the  program  is  the  scale  factors  for  the 
linear  combination  of  models.  What  each  models'  scale  factors 
represent  depends  on  the  normalization  dimensions  used  in 
creating  the  corresponding  synthesized  spectrum.  The  model 
spectra  do  not  all  have  the  same  dimensional  units.   In  the 
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case  of  atomic  spectra,  the  scale  factors  are  absolute 
intensity  measured  in  Rayleighs.  For  the  VK  and  LBH  bands  the 
scale  factors  are  a  measure  of  relative  intensity  and  must  be 
converted  to  absolute  values  to  obtain  an  intensity  profile. 
The  Nitric  Oxide  scale  factors  are  column  density,  in 
molecules  per  square  centimeter. 
Jb.  Synthetic  Spectra. 

Previous  analysis  of  the  March  1990  MUSTANG  data  includes: 
an  analysis  of  NO  at  2000-2500A  (Clayton,  1990)  ;  analysis  of  01 
2972A  and  Oil  2470A  emissions  (Andersen, 1990)  ;  and  a  study  of 
LBH  and  NO  at  1800-2100A  (Mack, 1991).  An  analysis  combining 
all  the  contributing  spectra  and  spanning  the  full  MUSTANG 
wavelength  range  had  not  been  performed,  but  previous  work 
provided  many  of  the  computer  programs  as  well  as  "initial 
guesses"  for  the  fitting  routine.  Other  models  and  programs 
arose  as  a  result  of  the  fitting  analysis.  Computer  models 
written  by  Cleary  (1986)  were  used  for  the  generation  of 
Nitrogen  Lyman-Birge-Hopf ield  and  Nitric  Oxide  7,  8,  and  e 
band  emissions.  The  Nitrogen  Vegard-Kaplan  model  was  written 
by  Siskind  and  Barth  (1987) . 

The  synthetic  spectra  for  individual  atomic  emissions  were 
created  by  convolving  appropriate  impulse  functions  with  the 
slit  function  of  the  instrument.  As  already  discussed,  in  the 
case  of  the  N+  2143A  doublet  a  pair  of  impulse  functions  in 
the  ratio  100:58  was  used.   The  previously  identified  atomic 
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contributions  to  the  MUSTANG  spectrum  include  01  2972 A,  Oil 
2470A,  and  Nil  2143A. 

One  other  atomic  spectrum,  a  prominent  line  emission  at 
2853 A  at  higher  altitudes,  has  been  included  in  the  fitting 
routine  as  a  result  of  observation.  Elements  with  emissions 
at  that  wavelength  include  nickel,  copper,  manganese,  iron, 
sodium,  cobalt,  and  oxygen.  Of  these,  oxygen  is  the  most 
abundant  in  the  atmosphere  where  the  emission  is  seen,  above 
200km.  Doubly  ionized  oxygen  has  an  emission  at  2853A  that 
although  not  previously  identified  in  dayglow  data,  is 
classified  "predicted"  by  Kelly  (1979) . 

A  collection  of  atomic  emissions  observed  during  the 
fitting  process,  but  not  verified  by  other  literature,  was 
classified  as  a  separate  spectrum.  As  discussed  in  Chapter 
II,  an  atomic  emission  is  a  single  line  vice  the  "hump"  of  a 
molecular  transition.  Even  after  convolution  with  the 
instrument  slit  function,  the  form  of  an  atomic  emission 
feature  is  distinctive  from  a  molecular  band.  There  were 
several  areas  in  the  data  which  matched  the  form  at  a 
sufficient  number  of  altitudes  as  to  suggest  an  atomic 
emission.  These  wavelengths  were  then  checked  against  tables 
for  possible  origins  and  included  in  the  synthetic  spectrum  if 
a  likely  emission  existed,  based  on  atmospheric  composition 
and  transition  probability.  The  resultant  synthetic  spectrum 
(designated  "Other  Atomics")  is  shown  in  Figure  5-1.  The 
program  used  to  create  this  spectrum,  along  with  possible 
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Figure  5-1   "Other  Atomics"  synthetic  spectrum 
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identifications  for  each  emission,  is  in  Appendix  B.  None  of 
these  emissions  occurred  with  sufficient  intensity  and 
regularity  to  merit  being  considered  as  a  separate  parameter. 
They  are  included  only  to  facilitate  the  fitting  procedure. 

2  .   General  Programming  Technique 
a.  Fitting  Routine 

The  GRIDLS  program  which  successfully  handled  72  points  of 
data  with  three  parameters  (Mack,  1991)  was  overloaded  by  the 
inclusion  of  all  512  points  of  wavelength  and  nine  parameters 
(i.e.,  the  scale  factors  for  V-K,  01,  Oil,  OIII,  LBH,  NO,  Nil, 
Other  Atomics,  and  Background) .  Several  modifications  were 
made  to  improve  the  data  handling  capabilities  of  the  program, 
and  other  characteristics  were  accepted  as  limitations. 

Previously,  the  program  would  allow  negative  scale 
factors,  effectively  allowing  subtraction  of  a  spectrum.  It 
also  would  enter  an  infinite  loop  if  the  scale  factor  equalled 
zero.  These  problems  were  eliminated  in  the  program  by 
requiring  that  zero  be  the  minimum  scale  factor  allowed,  and 
by  performing  a  simple  global  minimum  check  to  prevent  a  zero 
scale  factor  from  causing  an  infinite  loop. 

The  program  had  great  difficulty  converging  to  a  value  if 
the  "initial  guess"  for  the  scale  factor  was  far  from  the 
minimum  deviation  value.  This  problem  was  not  corrected,  but 
a  test-plotting  program  was  written  to  allow  one  to  visually 
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match  the  synthetic  curve  to  the  data  so  as  to  obtain 
reasonable  first  guess  values. 

The  standard  function  of  instrument  response  was  an  11- 
element  symmetric  slit  function  (Clayton,  1990;  Andersen, 
1990;  Mack,  1991) .  When  convolved  with  a  512-element 
synthetic  spectrum,  this  slit  function  would  sufficiently 
represent  the  instrument's  response  to  the  modeled  molecular 
spectrum.  However,  in  the  case  of  an  atomic  emission,  which 
occurs  as  a  line  at  a  particular  wavelength  vice  a  band  of 
lines,  the  slit  function  could  shift  the  apparent  peak  as 
much  as  1.7 A.  The  symmetric  function  centered  the  response  on 
the  applied  pixel,  a  3. 133 A  range,  instead  of  at  the 
particular  desired  wavelength.  A  new  set  of  five  slit 
functions,  shown  in  Figure  5-2,  was  developed  to  correctly 
model  atomic  emissions,  allowing  for  emissions  which  are  not 
centered  on  a  pixel. 

A  major  limitation  of  the  previous  version  of  the  GRIDLS 
fitting  routine  was  that  there  existed  no  absolute  measurement 
of  "goodness  of  fit."  Although  the  program  minimized  a  x2 
value,  the  values  were  not  weighted  and  thus  could  not  be 
compared  between  runs  or  altitudes  to  determine  goodness  of 
fit  as  parameters  were  tweaked.  Creating  a  proper  weighting 
for  the  data  evaluation  so  that  these  comparisons  could  be 
made  was  a  major  task. 
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b.    Weighting  Program 

The  GRIDLS  program  allows  for  statistical  weighting  of  the 
data,  instrumental  weighting,  or  no  weighting.  Instrumental 
weighting  was  chosen  for  its  advantage  of  stressing  datapoints 
which  have  small  uncertainties.  The  first  step  in  producing 
instrumental  data  weights  was  to  determine  the  standard 
deviations  of  the  mean  (sdom)  for  every  data  point  in  each 
altitude  bin.  A  second  source  of  instrumental  uncertainty  was 
instrument  response.  It  was  determined  that  the  detector  had 
a  linear  response  for  most  output  values,  but  reached 
saturation  above  900  decimal  units  (divided  by  instrument 
sensitivity  to  get  R/A) ,  so  data  points  in  the  saturation 
range  were  given  weights  based  on  the  sdom  being  egual  to  the 
data  value  itself.  The  exception  to  this  was  when  the  data 
value  was  so  small  as  to  give  a  large  weight  to  the  point 
(since  weight  oc  sdom"1)  ,  in  which  case  the  weight  was  assigned 
1/10  the  largest  weight  factor  for  that  altitude  (Cleary, 
private  communication,  1991)  .  Additionally,  the  data  points 
which  were  merely  averages  of  adjoining  values  due  to  lost 
data  points  were  given  weights  based  on  the  sdom  being  equal 
to  the  data  value.  Since  the  sdom  and  saturation  points  are 
different  for  each  altitude,  a  set  of  values  for  weighting  of 
the  512  data  points  was  stored  for  each  data  spectrum. 

The  data  pixels  corresponding  to  nitrogen  Second  Positive 
emissions  had  no  synthetic  input  being  matched  to  them.  The 
pixel  corresponding  to  each  wavelength  of  these  emissions  as 
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determined  by  Barth  (1965)  ,  plus  a  buffer  of  ±1  pixel  to 
account  for  instrument  response,  were  given  essentially  zero 
weight.  The  final  weight  program,  designated  WEIGHTS2,  is 
Appendix  C. 

B.   RESULTS 

1.   Comparison  of  Data  and  Fit 

a.  Nitric  Oxide  Bands 

In  the  plots  of  the  fits,  nitric  oxide  is  viewed  as  one 
entity  rather  than  plotting  the  y,  6,  and  e  bands  separately. 
This  is  for  ease  in  reading  the  plots.  Figure  5-3  shows  the 
components  of  the  nitric  oxide  synthetic  spectrum  for  115km. 
The  bands  are  not  equal  contributors,  as  the  S  band  has  a  25% 
fluorescent  efficiency  factor. 

Figure  5-4  is  a  plot  of  the  observed  data  and  its 
corresponding  fit.  The  data  is  indicated  by  the  solid  line, 
and  the  fit  obtained  from  the  computer  procedure  by  a  dotted 
line.  The  dashed  line  is  that  portion  of  the  synthetic  fit 
which  is  nitric  oxide.  From  the  plots  in  Figure  5-4,  one  can 
see  the  importance  of  nitric  oxide  with  decreasing  altitude. 
In  the  105km  plot,  nitric  oxide  accounts  for  almost  all  of  the 
dayglow  from  1800-2600A,  but  at  295km  the  NO  is  less 
prevalent. 

b.  Nitrogen  Bands 

Figure  5-5  shows  an  intensity  plot  (in  R/A)  of  the  LBH 
contribution  to  the  fit,  along  with  the  observed  spectrum,  for 
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Figure  5-3   Nitric  oxide  synthetic  bands  at  115km. 
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Figure  5-5   Comparison  of  data  &  LBH  synthetic  spectrum  for 
295,  155,  105km. 
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three  different  altitudes.  From  these  plots,  one  can  see  that 
the  Lyman-Birge-Hopf ield  band  of  N2  emissions  dominates  the 
wavelength  range  1800-2000A. 

Figure  5-6  shows  a  similar  altitude  progression  for  the 
Vegard-Kaplan  system,  emphasizing  the  (0,5)  transition  at 
2605A.  VK  emissions  are  prominent  from  2200-3300A,  and  in 
general  all  the  emission  features  show  the  same  characteristic 
growth  as  exhibited  by  the  (0,5)  transition. 
c.  Oxygen  Atomic  Emissions 

The  01  2972A  emission  as  fit  to  the  data  at  three 
different  altitudes  is  shown  in  Figure  5-7.  The  poor  fit  of 
this  feature  is  addressed  in  the  Discussion  section.  The  Oil 
2470A  line,  shown  in  Figure  5-8,  becomes  more  important  with 
increasing  altitude. 

2.   Profiles 

a.  Nitric  Oxide 

The  column  density  profile  for  Nitric  Oxide  is  shown  in 
Figure  5-9,  with  the  corresponding  values  given  in  Table  5-1. 
Values  obtained  from  upleg  data  are  indicated  with  squares  and 
downleg  values  by  triangles.  There  is  a  systematic 
uncertainty  of  15%.  The  values  obtained  show  general 
agreement  in  both  magnitude  and  slope  with  the  column 
densities  obtained  by  Mack  (1991) ,  Clayton  (1990)  and  Andersen 
(1990)  ,  all  of  whom  analyzed  subsets  of  the  MUSTANG  data.  The 
maximum  and  minimum  values  best  agreed  with  those  obtained  by 
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TABLE  5-1  NO  COLUMN  DENSITIES 


Mack  (1991) ,  with  the 
intermediate  values  best 
matching  the  fit  obtained  by 
Clayton  (1990) .  The 
significance  of  this  agreement 
is  that  a  full  spectrum 
fitting  procedure  (i.e., 
including  all  species  and  all 
wavelengths)  is  not  required 
to  accurately  determine  the 
Nitric  Oxide  density. 

b.    Nitrogen  Bands 

The  scale  factors  for  LBH 

bands    were    converted    to 

absolute  intensity  using  the 

algorithm    and    computer 


ALTITUDE(km 
) 

*downleg 

COLUMN  DENSITY 
(molecules/cm2) 

105* 

1.908e+16 

115* 

1.681e+16 

125* 

1.186e+16 

135* 

8.785e+15 

145* 

5.883e+15 

155 

2.427e+15 

155* 

3.596e+15 

165 

1.938e+15 

165* 

2.299e+15 

175 

1.779e+15 

175* 

1.520e+15 

185 

1.332e+15 

195 

9.850e+14 

205 

8.222e+14 

215 

6.321e+14 

225 

5.062e+14 

235 

4.003e+14 

245 

3.248e+14 

255 

2.554e+14 

265 

2.103e+14 

275 

1.665e+14 

285 

1.388e+14 

295 

1.158e+14 

305 

9.824e+13 

315 

6.240e+13 

programs  in  Appendix  D.  The  resultant  LBH  intensity  profile 
is  shown  in  Figure  5-10,  and  is  a  measurement  of  total  LBH 
intensity  as  it  includes  all  a->X  transitions.  The  intensity 
profile  for  VK  can  only  be  calculated  for  a  particular  v'->v" 
band,  and  Figure  5-11  shows  this  calculation  for  one  of  the 
more  prominent  systems,  the  (0,5). 
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c.    Atomic  Emissions 

The  intensity  profiles  for  01  2912k,     Oil  2470A,  OIII 
2853A,  and  Nil  2143A 

are  shown  in  Figure  TABLE  5-2  ATOMIC  INTENSITY  VALUES 
5-12,  with  the  values 
given  in  Table  5-2 . 
The  01  2972A  profile 
shows  a  distinctive 
peak  at  about  150km; 
the  0+  2470A  profile 
shows  a  broad  peak  at 
about  250km;  the  OIII 
profile  suggests  a 
peak  at  around  300km; 
and  the  Nil  reaches 
its  maximum  around 
150km. 

The  01  2972A  line 
emission  profile  was 


ALT(km) 
*down 
leg 

01  2972 A 

Oil  2470A 

OIII 
2853A 

Nil  2143A 
doublet 

105* 

1878.70 

7816.57 

0.00 

934.11 

115* 

1131.50 

6120.95 

0.00 

1113.95 

125* 

1207.23 

3477.89 

0.00 

1579.08 

135* 

1187.29 

1660.38 

0.00 

3106.81 

145* 

1652.85 

1633.97 

0.00 

2824.51 

155 

1777.33 

2065.98 

0.00 

1824.37 

155* 

1828.75 

1323.66 

0.00 

3469.86 

165 

1623.30 

1671.16 

0.00 

1715.49 

165* 

1692.26 

1264.87 

0.00 

2559.10 

175 

1351.22 

1436.91 

0.00 

1431.54 

175* 

1448.04 

1127.43 

0.00 

2175.71 

185 

1045.54 

1478.36 

0.00 

1243.45 

195 

907.09 

1535.49 

0.00 

1238.92 

205 

677.06 

1517.33 

0.00 

865.34 

215 

461.95 

1492.67 

0.00 

786.54 

225 

404.27 

1523.80 

0.00 

654.12 

235 

283.24 

1525.33 

0.00 

563.59 

245 

353.83 

1540.24 

0.00 

496.13 

255 

240.10 

1522.03 

0.00 

436.25 

265 

219.15 

1537.10 

22.56 

348.22 

275 

200.13 

1481.19 

25.02 

273.15 

285 

169.77 

1383.26 

25.48 

260.23 

295 

156.30 

1350.18 

30.51 

232.33 

305 

147.63 

1342.65 

32.31 

186.90 

315 

232.63 

1310.17 

29.27 

206.432 

found  to  have  very  little  agreement  with  previous  analysis. 
The  maximum  found  by  this  analysis  was  50%  of  that  found  when 
fitting  only  the  01  2972A  line  (Andersen,  1990) .  This  is  due 
to  the  uncertainty  of  the  wavelength  for  the  data,  since  two 
of  the  data  dropouts  were  in  this  region.   The  overall 
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uncertainty  in  the  01  profile  should  reflect  the  data  drops, 
thus  caution  is  urged  for  any  subsequent  analysis  of  this 
emission  profile. 
3.   Discussion 

a.  V-K  Quenching 

Figure  5-13  shows  an  overlay  of  the  nitrogen  LBH  and  VK 
profiles,  with  the  VK  scaled  up  by  a  factor  of  9.5  for 
comparison.  The  most  interesting  aspect  of  the  nitrogen 
profiles  is  that  they  have  roughly  the  same  slope  above  200 
kilometers,  but  the  VK  intensity  drops  off  rapidly  below  this 
height.  This  is  due  to  the  quenching  effect  of  N2,  NO,  0,  02, 
and  N  at  the  lower  altitudes.  From  auroral  studies,  atomic 
oxygen  is  thought  to  be  the  primary  quenching  agent  with 
molecular  oxygen  playing  a  role  below  135km  (Sharp,  1971) .  As 
previously  discussed,  the  upper  state  for  the  VK  emissions  is 
the  N2  A-state.  The  A-state  to  ground-state  energy  for  N2 
corresponds  to  an  excitation  energy  for  atomic  oxygen.  Below 
200  km,  the  O  density  is  high  enough  that  during  their 
metastable  lifetime  in  the  excited  state,  N2  molecules  collide 
with  and  excite  oxygen,  thus  losing  their  energy  without 
emitting  a  photon.  This  is  believed  to  be  the  first 
simultaneous  observation  of  these  two  profiles  in  the  dayglow. 

b.  Adjustments   to  Franck-Condon  Factors 

The  first  set  of  spectra  to  be  fitted  to  the  data  included 
Nitric  Oxide  y,    6,    and  e  bands,  Nitrogen  LBH  and  V-K  bands,  01 
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Figure  5-13   Intensity  plots  for  LBH  &  VK  showing  quenching 
of  VK. 
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2912k,  Oil  2470A,  and  OIII  2853A  emissions.  The  first  fine- 
tuning  of  the  fit  was  accomplished  by  plotting  the  fit  and 
data  simultaneously  in  4  00A  segments  and  looking  at  each  band 
for  trends  in  undershooting  or  overshooting  at  a  particular 
altitude,  or  for  all  altitudes  at  a  particular  wavelength. 
Any  trend  was  then  adjusted  up  or  down  in  the  original  model 
spectrum,  a  procedure  equivalent  to  varying  the  Franck-Condon 
factor  (the  factor  which  quantifies  probability  of  transition 
for  the  synthetic  models) .  This  process  was  very  tedious  and 
time  intensive,  but  a  second  fitting  yielded  an  average 
decrease  of  3  0%  in  x2 • 

After  the  first  set  of  changes  had  been  made  to  the 
Franck-Condon  factors,  the  same  process  of  plotting  and 
observing  trends  was  repeated.  Several  unidentified  peaks 
were  seen  in  the  data  which  are  believed  to  be  unidentified 
atomic  emissions,  since  their  profiles  matched  that  of  an 
impulse  function  convolved  with  the  instrument  slit  function. 
Other  than  the  well-documented  Nil  2143A  doublet,  which  was 
added  as  a  synthetic  spectrum,  these  emissions  had  not  been 
previously  reported  in  the  literature,  so  the  "Other  Atomics" 
spectrum  was  created  (see  discussion  above) .  The  two 
parameters,  Nil  and  "Other",  were  included  in  the  fitting 
routine.  Once  these  parameters  were  included,  the  adjustments 
to  the  Franck-Condon  factors  found  for  the  Nitric  Oxide  bands' 
synthetic  spectra  confirmed  those  used  by  Mack  (1991)  and 
Clayton  (1990) .   The  adjustments  to  the  Franck-Condon  factors 
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TABLE  5-3 


Adjustments  to  Franck-Condon  Factors 

Vegard-Kaplan 

Lyman-Birge- 
Hopfield 

Nitric  Oxide  y 

Nitric  Oxide  £ 

v',v" 

factor 

v',v" 

factor 

v',v" 

factor 

v'  ,v" 

factor 

0,k 
k=0...9 

2.6 

1,7 

1.7 

0,k 
k=0...8 

1.2 

0,0 

.63 

I.J 

j=0...9 

1.5 

3,9 

1.203 

1,0 

1.125 

0,1 

.9 

6,9 

.83 

5,11 

1.4 

1,1 

1.3 

0,2 

1.03 

4,8 

.67 

1.8 

1.3 

1.2 

1.8 

0,3 

1.11 

7,10 

.6 

3,10 

.8 

1,3 

1.6 

11,13 

.7 

4,11 

.9 

I.J 

j=4...10 

1.3 

12, 14 

.67 

5,12 

.757 

2,0 

.75 

2,8 

.67 

3,11 

.8 

2,2 

1.26 

5,10 

.5 

4,12 

.9 

2,3 

1.6 

8,12 

.5 

5,13 

1.2 

are  summarized  in  Table  5-3.  These  changes  to  the  accepted 
probabilities  of  transition  are  believed  to  be  actual 
corrections  and  not  just  a  function  of  the  instrumentation. 
It  is  expected  that  further  independent  measurements  will 
verify  these  corrections. 

The  intensity  profiles  were  derived  from  the  final 
parameters  of  the  fitting  routine,  but  analysis  of  the  data 
includes  looking  not  only  at  the  profiles  but  at  the  fits 
themselves.  A  complete  set  of  plots  of  the  data  showing  a 
comparison  of  data  and  fits,  as  well  as  accompanying  plots 
showing  the  breakdown  of  the  fits  by  constituent,  are  in 
Appendix  E.  Figures  5-14  and  5-15  show  sample  fits  of  the 
entire  wavelength  range  at  low  altitude  and  high  altitude.   It 
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Figure  5-14   Data  and  best-fitting  synthetic  spectrum  at 
13  5km. 
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Figure  5-15   Data  and  best-fitting  synthetic  spectrum  at 
305km. 
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is  apparent  from  these  figures  that  the  overall  fit  was 
excellent. 

It  should  be  pointed  out  that  there  is  a  disagreement 
between  the  data  and  fit  at  the  wavelengths  where  the  nitrogen 
Second  Positive  emissions  are  present.  Additionally,  the  only 
other  poorly  fit  feature  was  the  2930A  to  2990A  range  at 
altitudes  above  200  km.  This  is  one  of  the  areas  one  would 
like  to  be  able  to  fit  well,  since  it  includes  the  01  2912k 
emission.  Unfortunately,  two  of  the  dropped  data  points  are 
at  2969. ik  and  2972. 9A,  exactly  where  the  peak  of  01  is 
expected.  Also,  there  seems  to  be  another  emission  at  about 
2955A  which  is  not  explained  by  the  VK  bands  or  01  emission, 
creating  three  peaks  in  the  data  where  theory  predicts  only 
two.   This  feature  is  shown  in  Figure  5-16. 
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Figure  5-16  Averaged  spectrum  for  285km.  The  peak  at  2940A 
is  VK(0,7),  the  peak  at  2980A  is  due  to  01  2912k,  and  the 
feature  at  2955A  is  unidentified. 
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c.    Temperature  Profile 

In  addition  to  the  column  density  profile,  intensity 
profiles,  and  fits  of  model  spectra  to  data,  a  temperature 
profile  also  was  determined.  The  temperatures  described  by 
Mack  (1991)  were  used  as  initial  guesses.  Several  of  the 
temperatures  at  the  lowest  altitudes  were  increased  during 
fitting  to  improve  the  match  of  synthetic  spectra  to  data. 
The  resultant  temperature  profile  is  shown  in  Figure  5-17. 
The  temperatures  determined  by  the  fitting  procedure  are  all 
considerably  hotter  than  predicted  by  MSIS.  In  particular, 
the  asymptotic  temperature  is  400°  hotter.  This  is  possibly 
due  to  a  magnetic  storm  which  occurred  about  three  hours  prior 
to  launch. 
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Figure  5-17  The  temperatures  used  by  the  analysis  are 
indicated  by  symbols.  The  solid  line  is  the  MSIS 
temperature  profile  for  the  day  of  the  launch. 
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VI.   CONCLUSIONS 

A.   SUMMARY  OF  FINDINGS 

Thermospheric  dayglow  data  from  a  rocket-borne 
spectrograph  were  analyzed  in  this  study.  Computer 
synthesized  spectra  were  created  for  the  N2  Lyman-Birge- 
Hopfield  and  Vegard-Kaplan  bands;  NO  y,  6,  and  e  bands;  01 
2972A,  Oil  2470A,  and  OIII  2853A  atomic  lines;  and  Nil  2143A 
doublet  feature.  A  spectrum  of  other  possible  atomic 
emissions  and  a  background  factor  were  also  considered  in  the 
fit.  These  synthetic  models  were  fitted  to  the  observed  data, 
then  the  scaling  factors  were  used  to  determine  an  intensity 
profile  (column  density  profile  in  the  case  of  nitric  oxide) . 

The  intensity  and  column  density  profiles  obtained  were 
within  10-20%  of  the  values  acguired  by  simpler  fitting 
routines  in  all  cases  except  the  01  2972A  emission.  Thus  the 
analysis  supports  the  concept  of  fitting  of  just  one  feature 
to  arrive  at  density  and  intensity  profiles.  A  relationship 
was  noted  between  VK  and  LBH  emissions  that  would  seem  to 
support  the  theory  of  quenching  of  the  nitrogen  VK  system. 
The  data  above  roughly  3000A  could  not  be  modelled  without  a 
N2  Second  Positive  synthetic  spectrum.  The  loss  of  datapoints 
was  a  major  detriment  in  fitting  the  01  2972A  emission,  as 
there  were  two  data  dropouts  where  a  peak  was  expected.   The 
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temperature  profile  which  best  fit  the  data  was  higher  at  low 
altitudes  than  previous  analysis  had  predicted. 

B.   TOPICS  FOR  FURTHER  INVESTIGATION 

With  another  launch  scheduled  for  March  1992,  further 
spectral  analysis  of  the  March  1990  data  is  imprudent. 
Contemplating  the  prospective  data,  a  synthetic  spectrum  for 
the  N2  Second  Positive  would  enhance  the  analysis.  It  is 
expected  that  there  will  be  no  problem  with  data  dropouts  on 
the  next  flight  (Quint,  1991) ,  so  a  repeat  of  this  analysis 
(with  the  new  data)  will  greatly  improve  the  accuracy  of 
fitting  of  the  01  emission  and  possibly  show  that  it  can  be 
accurately  fitted  over  a  limited  wavelength  range. 

The  guenching  of  VK  below  200km  should  be  paralleled  by 
the  density  of  the  guenching  agent,  as  well  as  a  rise  in 
emissions  corresponding  to  the  relaxation  of  the  guenching 
agent  from  its  excited  state.  The  guenching  phenomenon  showed 
promise  for  rewarding  study,  but  was  not  the  subject  of  this 
thesis. 

Looking  toward  the  future,  the  MUSTANG  will  be  part  of  an 
operational  satellite  in  the  late  1990' s,  providing  a  constant 
flow  of  data  for  analysis.  Considering  the  numerous  new 
mathematical  theories  for  computer  network  calculations,  there 
may  be  a  method  better  suited  for  fitting  large  guantities  of 
data  than  the  GRIDLS  routine.  A  modification  to  the  fitting 
programs  to  place  less  stringent  reguirements  on  the  initial 
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programs  to  place  less  stringent  reguirements  on  the  initial 
guess  and  to  provide  a  more  formal  check  for  global  minimum 
would  greatly  improve  the  effectiveness  of  the  existing 
routine.  Additionally,  automation  of  the  analysis,  including 
derivation  of  the  electron  density  profile,  would  certainly 
add  to  the  viability  of  MUSTANG  as  a  provider  of  geophysical 
parameters  for  national  defense  agencies. 
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APPENDIX  A 

The  following  table  shows  the  unshifted  and  shifted  array 
elements  of  data,  and  their  corresponding  wavelengths. 


Record:  [mustang] m90  avg/um°0  avg 

Entry  M 
(pixel) 
Shifted 

Entry  # 

(pixel) 

Unshifted 

Wavelength 

A 

0 

0 

1,798.00 

1 

1 

1,801.13 

2 

2 

1,804.27 

3 

3 

1,807.40 

4 

4 

1,810.53 

5 

5 

1,813.67 

6 

6 

1,816.80 

7 

7 

1,819.93 

8 

8 

1,823.06 

9 

9 

1,826.20 

10 

10 

1,829.33 

11 

11 

1,832.46 

12 

12 

1,835.60 

13 

13 

1,838.73 

H 

14 

1,841.86 

15 

15 

1,845.00 

16 

1,848.13 

17 

16 

1,851.26 

18 

17 

1,854.39 

19 

18 

1,857.53 

20 

19 

1,860.66 

21 

20 

1,863.79 

22 

21 

1,866.93 

23 

22 

1,870.06 

24 

23 

1,873.19 

25 

24 

1,876.33 

26 

25 

1,879.46 

27 

26 

1,882.59 
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28 

27 

1,885.72 

29 

28 

1,888.86 

30 

29 

1,891.99 

31 

30 

1,895.12 

32 

31 

1,898.26 

33 

1,901.39 

34 

32 

1,904.52 

35 

33 

1,907.66 

36 

34 

1,910.79 

37 

35 

1,913.92 

38 

36 

1,917.05 

39 

37 

1,920.19 

40 

38 

1,923.32 

41 

39 

1,926.45 

42 

40 

1,929.59 

43 

41 

1,932.72 

44 

42 

1,935.85 

45 

43 

1,938.99 

46 

44 

1,942.12 

47 

45 

1,945.25 

48 

46 

1,948.38 

49 

47 

1,951.52 

50 

1,954.65 

51 

48 

1,957.78 

52 

49 

1,960.92 

53 

50 

1,964.05 

54 

51 

1,967.18 

55 

52 

1,970.32 

56 

53 

1,973.45 

57 

54 

1,976.58 

58 

55 

1,979.71 

59 

56 

1,982.85 

60 

57 

1,985.98 

61 

58 

1,989.11 

62 

59 

1,992.25 

63 

60 

1,995.38 

64 

61 

1,998.51 

59 


65 

62 

2,001.65 

66 

63 

2,004.78 

67 

2,007.91 

68 

64 

2,011.04 

69 

65 

2,014.18 

70 

66 

2,017.31 

71 

67 

2,020.44 

72 

68 

2,023.58 

73 

69 

2,026.71 

74 

70 

2,029.84 

75 

71 

2,032.98 

76 

72 

2,036.11 

77 

73 

2,039.24 

78 

74 

2,042.37 

79 

75 

2,045.51 

80 

76 

2,048.64 

81 

77 

2,051.77 

82 

78 

2,054.91 

83 

79 

2,058.04 

84 

2,061.17 

85 

2,064.31 

86 

80 

2,067.44 

87 

81 

2,070.57 

88 

82 

2,073.70 

89 

83 

2,076.84 

90 

84 

2,079.97 

91 

85 

2,083.10 

92 

86 

2,086.24 

93 

87 

2,089.37 

94 

88 

2,092.50 

95 

89 

2,095.64 

96 

90 

2,098.77 

97 

91 

2,101.90 

98 

92 

2,105.03 

99 

93 

2,108.17 

100 

94 

2,111.30 

101 

2,114.43 
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102 

2,117.57 

103 

95 

2,120.70 

104 

96 

2,123.83 

105 

97 

2,126.97 

106 

98 

2,130.10 

107 

99 

2,133.23 

108 

100 

2,136.36 

109 

101 

2,139.50 

110 

102 

2,142.63 

111 

103 

2,145.76 

112 

104 

2,148.90 

113 

105 

2,152.03 

114 

106 

2,155.16 

115 

107 

2,158.30 

116 

108 

2,161.43 

117 

109 

2,164.56 

118 

110 

2,167.69 

119 

2,170.83 

120 

111 

2,173.96 

121 

112 

2,177.09 

122 

113 

2,180.23 

123 

114 

2,183.36 

124 

115 

2,186.49 

125 

116 

2,189.63 

126 

117 

2,192.758 

127 

118 

2,195.891 

128 

119 

2,199.024 

129 

120 

2,202.157 

130 

121 

2,205.290 

131 

122 

2,208.423 

132 

123 

2,211.556 

133 

124 

2,214.689 

134 

125 

2,217.822 

135 

126 

2,220.955 

136 

2,224.088 
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137 

127 

2,227.221 

138 

128 

2,230.354 

139 

129 

2,233.487 

140 

130 

2,236.620 

141 

131 

2,239.753 

142 

132 

2,242.886 

143 

133 

2,246.019 

144 

134 

2,249.152 

145 

135 

2,252.285 

146 

136 

2,255.418 

147 

137 

2,258.551 

148 

138 

2,261.684 

149 

139 

2,264.817 

150 

140 

2,267.950 

151 

141 

2,271.083 

152 

142 

2,274.216 

153 

2,277.349 

154 

143 

2,280.482 

155 

144 

2,283.615 

156 

145 

2,286.748 

157 

146 

2,289.881 

158 

147 

2,293.014 

159 

148 

2,296.147 

160 

149 

2,299.280 

161 

150 

2,302.413 

162 

151 

2,305.546 

163 

152 

2,308.679 

164 

153 

2,311.812 

165 

154 

2,314.945 

166 

155 

2,318.078 

167 

156 

2,321.211 

168 

157 

2,324.344 

169 

158 

2,327.477 

170 

2,330.610 

171 

159 

2,333.743 

172 

160 

2,336.876 

173 

161 

2,340.009 
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174 

162 

2,343.142 

175 

163 

2,346.275 

176 

164 

2,349.408 

177 

165 

2,352.541 

178 

166 

2,355.674 

179 

167 

2,358.807 

180 

168 

2,361.940 

181 

169 

2,365.073 

182 

170 

2,368.206 

183 

171 

2,371.339 

184 

172 

2,374.472 

185 

173 

2,377.605 

186 

174 

2,380.738 

187 

2,383.871 

188 

175 

2,387.004 

189 

176 

2,390.137 

190 

177 

2,393.270 

191 

178 

2,396.403 

192 

179 

2,399.536 

193 

180 

2,402.669 

194 

181 

2,405.802 

195 

182 

2,408.935 

196 

183 

2,412.068 

197 

184 

2,415.201 

198 

185 

2,418.334 

199 

186 

2,421.467 

200 

187 

2,424.600 

201 

188 

2,427.733 

202 

189 

2,430.866 

203 

190 

2,433.999 

204 

2,437.132 

205 

191 

2,440.265 

206 

192 

2,443.398 

207 

193 

2,446.531 

208 

194 

2,449.664 

209 

195 

2,452.797 

210 

196 

2,455.930 
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211 

197 

2,459.063 

212 

198 

2,462.196 

213 

199 

2,465.329 

214 

200 

2,468.462 

215 

201 

2,471.595 

216 

202 

2,474.728 

217 

203 

2,477.861 

218 

204 

2,480.994 

219 

205 

2,484.127 

220 

206 

2,487.260 

221 

2,490.393 

222 

207 

2,493.526 

223 

208 

2,496.659 

224 

209 

2,499.792 

225 

210 

2,502.925 

226 

211 

2,506.058 

227 

212 

2,509.191 

228 

213 

2,512.324 

229 

214 

2,515.457 

230 

215 

2,518.590 

231 

216 

2,521.723 

232 

217 

2,524.856 

233 

218 

2,527.989 

234 

219 

2,531.122 

235 

220 

2,534.255 

236 

221 

2,537.388 

237 

222 

2,540.521 

238 

2,543.654 

239 

2,546.787 

240 

223 

2,549.920 

241 

224 

2,553.053 

242 

225 

2,556.186 

243 

226 

2,559.319 

244 

227 

2,562.452 

245 

228 

2,565.585 

246 

229 

2,568.718 

247 

230 

2,571.851 
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248 

231 

2,574.984 

249 

232 

2,578.117 

250 

233 

2,581.250 

251 

234 

2,584.383 

252 

235 

2,587.516 

253 

236 

2,590.649 

254 

237 

2,593.782 

255 

238 

2,596.915 

256 

2,600.048 

257 

239 

2,603.181 

258 

240 

2,606.314 

259 

241 

2,609.447 

260 

242 

2,612.580 

261 

243 

2,615.713 

262 

244 

2,618.846 

263 

245 

2,621.979 

264 

246 

2,625.112 

265 

247 

2,628.245 

266 

248 

2,631.378 

267 

249 

2,634.511 

268 

250 

2,637.644 

269 

251 

2,640.777 

270 

252 

2,643.910 

271 

253 

2,647.043 

272 

2,650.176 

273 

254 

2,653.309 

274 

255 

2,656.442 

275 

256 

2,659.575 

276 

257 

2,662.708 

277 

258 

2,665.841 

278 

259 

2,668.974 

279 

260 

2,672.107 

280 

261 

2,675.240 

281 

262 

2,678.373 

282 

263 

2,681.506 

283 

264 

2,684.639 
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284 

265 

2,687.772 

285 

266 

2,690.905 

286 

267 

2,694.038 

287 

268 

2,697.171 

288 

269 

2,700.304 

289 

2,703.437 

290 

270 

2,706.570 

291 

271 

2,709.703 

292 

272 

2,712.836 

293 

273 

2,715.969 

294 

274 

2,719.102 

295 

275 

2,722.235 

296 

276 

2,725.368 

297 

277 

2,728.501 

298 

278 

2,731.634 

299 

279 

2,734.767 

300 

280 

2,737.900 

301 

281 

2,741.033 

302 

282 

2,744.166 

303 

283 

2,747.299 

304 

284 

2,750.432 

305 

285 

2,753.565 

306 

2,756.698 

307 

286 

2,759.831 

308 

287 

2,762.964 

309 

288 

2,766.097 

310 

289 

2,769.230 

311 

290 

2,772.363 

312 

291 

2,775.496 

313 

292 

2,778.629 

314 

293 

2,781.762 

315 

294 

2,784.895 

316 

295 

2,788.028 

317 

296 

2,791.161 

318 

297 

2,794.294 

319 

298 

2,797.427 

320 

299 

2,800.560 
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321 

300 

2,803.693 

322 

301 

2,806.826 

323 

2,809.959 

324 

302 

2,813.092 

325 

303 

2,816.225 

326 

304 

2,819.358 

327 

305 

2,822.491 

328 

306 

2,825.624 

329 

307 

2,828.757 

330 

308 

2,831.890 

331 

309 

2,835.023 

332 

310 

2,838.156 

333 

311 

2,841.289 

334 

312 

2,844.422 

335 

313 

2,847.555 

336 

314 

2,850.688 

337 

315 

2,853.821 

338 

316 

2,856.954 

339 

317 

2,860.087 

340 

2,863.220 

341 

318 

2,866.353 

342 

319 

2,869.486 

343 

320 

2,872.619 

344 

321 

2,875.752 

345 

322 

2,878.885 

346 

323 

2,882.018 

347 

324 

2,885.151 

348 

325 

2,888.284 

349 

326 

2,891.417 

350 

327 

2,894.550 

351 

328 

2,897.683 

352 

329 

2,900.816 

353 

330 

2,903.949 

354 

331 

2,907.082 

355 

332 

2,910.215 

356 

333 

2,913.348 

357 

2,916.481 
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.358 

334 

2,919.614 

359 

335 

2,922.747 

360 

336 

2,925.880 

361 

337 

2,929.013 

362 

338 

2,932.146 

363 

339 

2,935.279 

364 

340 

2,938.412 

365 

341 

2,941.545 

366 

342 

2,944.678 

367 

343 

2,947.811 

368 

344 

2,950.944 

369 

345 

2,954.077 

370 

346 

2,957.210 

371 

347 

2,960.343 

372 

348 

2,963.476 

373 

349 

2,966.609 

374 

2,969.742 

375 

2,972.875 

376 

350 

2,976.008 

377 

351 

2,979.141 

378 

352 

2,982.274 

379 

353 

2,985.407 

380 

354 

2,988.540 

381 

355 

2,991.673 

382 

356 

2,994.806 

383 

357 

2,997.939 

384 

358 

3,001.072 

385 

359 

3,004.205 

386 

360 

3,007.338 

387 

361 

3,010.471 

388 

362 

3,013.604 

389 

363 

3,016.737 

390 

364 

3,019.870 

391 

365 

3,023.003 
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392 

3,026.136 

393 

366 

3,029.269 

394 

367 

3,032.402 

395 

368 

3,035.535 

396 

369 

3,038.668 

397 

370 

3,041.801 

398 

371 

3,044.934 

399 

372 

3,048.067 

400 

373 

3,051.200 

401 

374 

3,054.333 

402 

375 

3,057.466 

403 

376 

3,060.599 

404 

377 

3,063.732 

405 

378 

3,066.865 

406 

379 

3,069.998 

407 

380 

3,073.131 

408 

381 

3,076.264 

409 

3,079.397 

410 

382 

3,082.530 

411 

383 

3,085.663 

412 

384 

3,088.796 

413 

385 

3,091.929 

414 

386 

3,095.062 

415 

387 

3,098.195 

416 

388 

3,101.328 

417 

389 

3,104.461 

418 

390 

3,107.594 

419 

391 

3,110.727 

420 

392 

3,113.860 

421 

393 

3,116.993 

422 

394 

3,120.126 

423 

395 

3,123.259 

424 

396 

3,126.392 

425 

397 

3,129.525 

426 

3,132.658 

427 

398 

3,135.791 

428 

399 

3,138.924 
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429 

400 

3,142.057 

430 

401 

3,145.190 

431 

402 

3,148.323 

432 

403 

3,151.456 

433 

404 

3,154.589 

434 

405 

3,157.722 

435 

406 

3,160.855 

436 

407 

3,163.988 

437 

408 

3,167.121 

438 

409 

3,170.254 

439 

410 

3,173.387 

440 

411 

3,176.520 

441 

412 

3,179.653 

442 

413 

3,182.786 

443 

3,185.919 

444 

414 

3,189.052 

445 

415 

3,192.185 

446 

416 

3,195.318 

447 

417 

3,198.451 

448 

418 

3,201.584 

449 

419 

3,204.717 

450 

420 

3,207.850 

451 

421 

3,210.983 

452 

422 

3,214.116 

453 

423 

3,217.249 

454 

424 

3,220.382 

455 

425 

3,223.515 

456 

426 

3,226.648 

457 

427 

3,229.781 

458 

428 

3,232.914 

459 

429 

3,236.047 

460 

3,239.180 

461 

430 

3,242.313 

462 

431 

3,245.446 

463 

432 

3,248.579 

464 

433 

3,251.712 

465 

434 

3,254.845 
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466 

435 

3,257.978 

467 

436 

3,261.111 

468 

437 

3,264.244 

469 

438 

3,267.377 

470 

439 

3,270.510 

471 

440 

3,273.643 

472 

441 

3,276.776 

473 

442 

3,279.909 

474 

443 

3,283.042 

475 

444 

3,286.175 

476 

445 

3,289.308 

477 

3,292.441 

478 

446 

3,295.574 

479 

447 

3,298.707 

480 

448 

3,301.840 

481 

449 

3,304.973 

482 

450 

3,308.106 

483 

451 

3,311.239 

484 

452 

3,314.372 

485 

453 

3,317.505 

486 

454 

3,320.638 

487 

455 

3,323.771 

488 

456 

3,326.904 

489 

457 

3,330.037 

490 

458 

3,333.170 

491 

459 

3,336.303 

492 

460 

3,339.436 

493 

461 

3,342.569 

494 

3,345.702 

495 

462 

3,348.835 

496 

463 

3,351.968 

497 

464 

3,355.101 

498 

465 

3,358.234 

499 

466 

3,361.367 

500 

467 

3,364.500 

501 

468 

3,367.633 

502 

469 

3,370.766 
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503 

470 

3,373.899 

504 

471 

3,377.032 

505 

472 

3,380.165 

506 

473 

3,383.298 

507 

474 

3,386.431 

508 

475 

3,389.564 

509 

476 

3,392.697 

510 

477 

3,395.830 

511 

3,398.963 
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APPENDIX  B 

Appendix  B  contains  the  least  squares  fitting  program  and 
its  subroutines.  Additionally,  the  synthetic  spectrum 
generation  routines  are  given. 
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;Walden,1991,  modified  from  Mack, 1991 

@getf unc 

@f chisqr 

pro  nimpfit_all 

This  program  fits  data  for  the  1800  to  3400A  range  with  models  of  the 
NO  gamma,  delta,  and  epsilon  bands,  N2  LBH  emissions,  N2  Vegard- 
Kaplan  bands,  oxygen  emissions  at  2470,  2972,  &  2853A,  N+  doublet 
at  2143A,  various  atomic  emissions,  and  the  background  level 
through  an  iterative  process.  It  uses  gridls.pro,  getfunc.pro,  and 
fchisqr.pro  as  subprograms;  it  uses  parameter.dat  as  an  in/out  file; 
and  it  writes  the  fit  by  altitude  to  nimpfit_all.dat.   The  subprograms 
utilize  synthetic  spectra  which  are  stored  in  allvk.dat,  all_lbh.dat, 
nii_2143doub.dat,  nodlt.dat,  nogammack.dat,  noeps.dat,  oiii_"2~853  .dat 
oii_2470.dat,  oi_2972.dat,  and  otherspec.dat. 

parameter.dat  -  file  of  parameters-  records  have  25  elements  each 

record!-  contents 

0  -  temperature 

1  -  altitude  record  number  (1  through  25;1-17  up, 18-25  down) 

2  -  chi  squared  of  the  fit 

3  -  scale  factor  for  0  2972A  emission,  in  Rayleighs ( intensity ) 

4  -  sigma  for  0  scale  factor 

5  -  scale  factor  for  0+  2470A  emission,  in  Rayleighs 

6  -  sigma  for  0+  scale  factor 

7  -  scale  factor  for  0++  2853A  emission,  in  Rayleighs 

8  -  sigma  for  0++  scale  factor 

9  -  scale  factor  for  the  N2  LBH  bands 

10  -  sigma  for  LBH  scale  factor 

11  -  scale  factor  for  the  NO  gamma,  delta,  and  epsilon  bands,  in 
molecules  per  square  centimeter  (column  density) 

12  -  sigma  for  NO  scale  factor 

13  -  scale  factor  for  the  N2  V-K  bands 

14  -  sigma  for  VK  scale  factor 

15  -  scale  factor  for  N+  2143A  emission,  in  Rayleighs 

16  -  sigma  for  N+  doublet  scale  factor 

17  -  scale  factor  for  synthetic  atomics  spectrum,  in  Rayleighs 

18  -  sigma  for  atomic  spectrum 

19  -  scale  factor  for  the  background,  in  Rayleighs  per  Angstrom 

20  -  sigma  for  background  scale  factor 

21  -  altitudes  in  kilometers 

para  -  array (( 2*nterms )+2 )  in  same  order  as  record  order  in  file 

a    -  array ( nterms+1 )  of  temp  followed  by  scale  factors 

yfit  -  512  element  array  of  best  fit 

start-  first  record  number 

stop  -  last  record  number 

openr , 30 ,' [mustang ]m90_int .dat '      ;open  data  file 
aa=assoc(30,fltarr(512) ) 

openu, 31 , ' nimpf i t_all . dat '  ;open  file  to  save  yfit 

bb=assoc(31,fltarr(512) ) 

openu, 32 ,' parameter .dat '  ;open  parameter  file 

pp=assoc( 32,fltarr(25) ) 

openr ,33, ' [walden.datafix ]m9  0_domwt .dat' 
dd=assoc( 33,fltarr(512)) 
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start=0 

stop=511 

npts=( stop-start )+l 

nterms  =  9 

para=fltarr(2*nterms+3) 
a=fltarr(nterms+l) 


#  of  points  to  be  fit 

#  of  terms  to  consider,  change  if 
parameters  are  added 


;fill  para  with  values  for  alt  i 


;assign  first  element  in  a  temp 
;fill  a  with  scale  factors 


for  i=0,24  do  begin 

for  j=0 , ( 2*nterms ) +2  do  begin 
hold=pp( j ) 
para( j )=hold( i ) 
endf or 
alt=para ( 1 ) 
count=0 
a( 0)=para( 0) 

for  j=3 , ( 2*nterms )+l ,+2  do  begin 
count=count+l 
a ( count )=para ( j ) 
endf or 
y=aa( start:stop,alt) 

;deltaa  is  step  size  for  estimations,  cannot  be  zero,  first  element  is  dummy 
h=a 
for  j=l, nterms  do  begin 

h( j)=a( j) 

if  (h(j)  eq  0.0)  then  h ( j ) =1000 . 00 
endf or 
deltaa=h/10. 


mode=l 
weights=dd( alt 


;use  sdom  for  weighting 


gridls,y,weights,npts , nterms ,mode , a ,deltaa , sigmaa ,yfit,chisqr, start, stop 
bb(alt)=yfit 

count=0  ;save  parameters  from  best  fit 

para( 2 )=chisqr 

for  j=3 , ( 2*nterms )+l ,+2  do  begin 

count=count+l 

para( j  )=a( count ) 

para( j+1 )=sigmaa( count ) 
endf or 
for  j=0 , ( 2*nterms )+2  do  begin 

hold=pp( j ) 

hold( i )=para ( j ) 

pp( j )=hold 
endf or 


;put  parameters  in  appropriate 
jrecord  in  parameter.dat 


ret 
end 
@gr 


urn 
idls 


print, 'alt=' ,alt, '  a= ' ,a, '  chisqr= ' ,chisqr 
endf or 

close , 30 
close , 31 
close , 32 
close , 33 
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;Walden,1991  modified  from  Mack, 1991 
function  getfunc,a 

This  function  computes  the  synthetic 
for  given  scaling  factors. 

INPUTS 
a    9-D  array  of  parameters 

-  temperature  at  altitude 

-  scale  factor  for 

-  scale  factor  for 

-  scale 

-  scale 

-  scale 


spectrum  from  1800  to  3400  A 


9-D 

a(0) 

a(l) 

a(2) 

a(3) 

a(4) 

a(5) 

a(6) 

a(7) 

a(8) 

a(9) 

factor 
factor 
factor 


of  interest 
0  emission  @  2972 
0+  emission  @  2470 
0++  emission  @  2853 
N2  LBH  bands 
the  sum  of  NO  gamma 


for 

for 

for  the  sum  of  NO  gamma,  delta, 
and  epsilon  bands 

scale  factor  for  N2  Vegard-Kaplan  bands 
scale  factor  for  N2  doublet  emission  @2143 
scale  factor  for  "Other  Atomics"  spectrum 
scale  factor  for  background 

GETfunc  returns  an  array  of  511  points  to  be  compared  to  the 
experimentally  recorded  profile  from  1800  to  3400  A. 


LBH 
NOSYN 

VK 

01 

Oil 

OIII 

N2 

AT 

openr , 20 
rr=assoc 
LBH=rr(a 

openr , 21 
ss=assoc 
NOGAM=ss 

openr , 22 
tt=assoc 
NODLT=tt 

openr , 23 
uu=assoc 
NOEPS=uu 

openr , 24 
vv=assoc 


The  N2  LBH  bands  convolved  with  the  slit  function 

The  sum  of  the  NO  gamma,  delta,  and  epsilon  bands  convolved 

with  the  slit  function 

The  N2  V-K  bands  convolved  with  the  slit  function 

The  01  emission  @  2972  A  convolved  with  the  slit  function 

The  Oil  emission  @  2470  A  convolved  with  the  slit  function 

The  OIII  emission  @  2853  A  convolved  with  the  slit  function 

The  N2  emission  (§2143  A  convolved  with  the  slit  function 

The  spectrum  of  possible  atomic  emissions 


'ALL_LBH.DAT' 
20,fltarr(512)  ) 
0)/100) 

'NOGAMMACK.DAT' 
21,fltarr(512)  ) 
a(0)/100) 

'NODLT.DAT' 
22,fltarr(512)  ) 
a(0)/100) 

'NOEPS.DAT' 
23,fltarr(512)  ) 
a(0)/100) 

'ALLVK.DAT' 
24,fltarr(512)  ) 


VK=vv(a(0)/100) 

openr, 25, 'OI_2972.DAT' 
ww=  assoc(25,fltarr(512) ) 
OI=ww( 0 ) 

openr, 26, 'Oil  2470.DAT' 
xx=assoc(26, fTtar r ( 512 ) ) 
OII-xx(O) 

openr, 27, 'OIII_2853.DAT' 
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yy=assoc( 27,fltarr(512) ) 
OIII=yy( 0) 

openr,28, 'NII_2143DOUB.DAT' 
zz=assoc(28,fltarr(512) ) 
N2=zz( 1 ) 

openr,29, 'otherspec.dat' 
atsp=assoc(29,fltarr(512) ) 
AT=atsp( 1 ) 

NOSYN  =   NOGAM  +  .25*NODLT  +  NOEPS  ;25%  flourescent  efficiency  of 

;delta  band 

getfunc  =  a(4)*LBH  +  a ( 5 ) *NOSYN/3 . 133e6  +  a(6)*VK  +  a(l)*OI  $ 

+  a(2)*OII  +  a(3)*OIII  +  a(7)*N2  +a(8)*AT  +  a(9)/3.133 
;background  is  divided  by  wavelength  bin  so  that 
;scale  factor  is  absolute  intensity,  NOSYN  is  divided 
;by  3.133e6  for  same  reason 

close , 20 
close , 21 
close , 22 
close , 23 
close , 24 
close , 25 
close , 26 
close , 27 
close , 28 
close , 29 

return, getfunc 
end 
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;Mack,1991 

function  fchisqr,y,weights,npts,free ,mode ,yfit, start, stop 
This  function  evaluates  the  reduced  chi  squared 
for  a  fit  to  the  data. 
Adapted  from  Bevington,  pg  194. 

INPUTS 
y      -  array  of  data  points 

weights-  array  of  weightings  based  on  sdom's  of  data  points 
npts    -  number  of  data  points  being  fit 
free    -  number  of  degrees  of  freedom 
mode    -  determines  method  of  weighting  the  fit 
+1  instrumental  weighting( sdom  weights) 

0  no  weighting 
-1  statistical  weighting 
yfit    -  array  of  fitted  points 
start   -  first  record  element  to  be  fit 
stop   -  last  record  element  to  be  fit 

count=0 
chisqr=0 . 

if  mode  eq  0  then  wt=l. 

if  mode  gt  0  then  wt=weights 

if  mode  It  0  then  wt=l./abs(y) 

for  i=start,stop  do  begin 

chisqr=chisqr  +  wt ( i )*( y( count )-yf it ( count )) "2 

count=count  +  1 
endf or 

fchisqr=chisqr/free 
return , f chisqr 
md 
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;Walden,1991 

oro  lbhcalc , wl , int , sli t 

openu,10, 'all_lbh.dat' 
lbh=assoc(10,fltarr(512) ) 


Calculates  synthetic  lbh  spectra  at  temps 
from  100  to  2000  deg  for  wavelengths  from 
1800  to  3400  A.  Results  are  convolved  with 
MUSTANG  instrument  slit  function. 


wl=indgen(512)*3.133  +  1798 


openr ,11, ' [walden]slit.daf 
sl=assoc(ll,fltarr(ll) ) 
slit-sl(l) 


i  =  0 
for 


=100,20 
=  i  +  l 
nt=wl*0 
ntl=wl* 
nt2=wl* 
nt3=wl* 
nt4=wl* 
nt5=wl* 
nt6=wl* 
nt7=wl* 
nt8=wl* 
nt9=wl* 
ntl0=wl 
bhsyn, t 
bhsyn, t 
bhsyn, t 
bhsyn, t 
nt=int+ 
bhsyn, t 
bhsyn, t 
bhsyn, t 
nt=int+ 
bhsyn, t 
bhsyn, t 
bhsyn, t 
nt=int+ 
bhsyn, t 
nt=int+ 
bhsyn, t 
bhsyn, t 
bhsyn, t 
nt=int+ 
bhsyn, t 
nt=int+ 
bhsyn, t 
bhsyn, t 
nt=int+ 
bhsyn, t 
bhsyn, t 
bhsyn, t 
bhsyn, t 
nt=int+ 
bhsyn, t 
nt=int+ 
bhsyn, t 
bhsyn, t 


00,100  do  begin 


Wavelength 

variables  for  tweaking  individual  transitions 


0 

4,9 

0,6 

5,1 

1,7 

ntl 

6,1 

2,8 

3,9 

ntl 

4,1 

0,7 

5,1 

nt2 

1,8 

nt3 

6,1 

2,9 

3,1 

nt4 

4,1 

nt5 

0,8 

5,1 

nt6 

1,9 

6,1 

2,1 

3,1 

nt7 

4,1 

nt8 

0,9 

5,1 


,wl, 

,wl, 

0,wl 

,wl, 

0*1. 

l,wl 

,wl, 

,wl, 

*1.2 

0,wl 

,wl, 

l,wl 

*1.4 

,wl, 

*1.3 

2,wl 

,wl, 

0,wl 

*.8 

l,wl 

*.9 

,wl, 

2,wl 

*.75 

,wl, 

3,wl 

0,wl 

l,wl 

*.8 

2,wl 

*.9 

,wl , 

3,wl 


int 

int 

,  int 

intlO 

7 

,  int 

int 

intl 

03 

,  int 

int 

,  int2 

int3 

,  int 
int 

,int4 

,  int5 

int 
,  int6 
7 
int 

,  int 
,  int 
,int7 

,  int8 

int 
,int9 


1801.0 
1804.6 
1817.9 
1820.8 

1835.0 
1837.2 
1853.8 

1870.8 
1877.7 
1887.9 

1894.2 


1905 
1910 
1927 


;  1945.0 

;  19  5  5.9 

;  1962.4 

;  1972.6 

;  1980.1 

;  1989.6 

;  2  0  0  6.8 

;  2024.2 

;  2  0  4  0 

;  2  0  41.9 
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int=int+ 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 

lbhsyn 


int9*l .2 


,t 

,1 

,10 

,wl, 

,t 

,6 

rl4 

rWl, 

ft 

r2 

,11 

,Wl, 

,t 

,3 

,12 

,wl, 

,t 

,4 

,13 

,  wl , 

,  t 

,5 

,14 

,wl, 

,t 

,0 

,10 

,wl, 

,t 

,6 

,15 

,wl, 

ft 

,1 

,11 

,  wl , 

,  t 

,2 

,12 

,wl, 

,  t 

,3 

,13 

,wl, 

,t 

,4 

,14 

,wl, 

,t 

,5 

,15 

,wl, 

ft 

r0 

,11 

,wl, 

,t 

,6 

,16 

,wl, 

,t 

,1 

,12 

,wl, 

ft 

,2 

,13 

rWl, 

,t 

,3 

,14 

,wl, 

,t 

,4 

,15 

,wl, 

,t 

,5, 

,16 

,wl, 

,  t 

-6, 

17 

,wl, 

,  t 

■  1, 

13, 

wl , 

,  t 

2, 

14, 

wl, 

ft, 

3, 

15, 

wl, 

ft, 

4, 

16, 

wl, 

ft, 

5, 

17, 

wl, 

ft, 

6, 

18, 

wl, 

ft, 

2, 

15, 

wl, 

ft, 

3, 

16, 

wl, 

ft, 

4, 

17, 

wl , 

ft, 

5, 

18, 

wl , 

ft, 

6, 

19< 

wl, 

ft, 

4, 

18, 

wl, 

ft, 

5, 

19, 

wl , 

ft, 

6, 

20, 

wl, 

ft, 

6, 

21, 

wl, 

int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 
int 


2056 
2059 
2073 
2091 
2108 
2126 
2130 
2144 
2147 
2164 
2181 
2198 
2217 
2227 
2235 
2244 
2261 
2279 
2297 
2314 
2332 
2349 
2366 
2384 
2401 
2419 
2437 
2479 
2496 
2514 
2531 
2549 
2636 
2653 
2670 
2801 


print, total ( int ) 
spec=convol(int,slit) 
lbh( i )=spec/total ( spec ) 
endf or 

close , 10 
close , 11 

eturn 

nd 

[SYNSPEC] LBHSYN 
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ilden,1991 

allvkcalc,wl, int, slit 
mr ,1, ' [WALDEN]slit .dat' 
■nu, 3 , ' allvk .dat' 
=assoc(l,fltarr(ll) ) 
t-sl(l) 

=assoc(3,fltarr(512) ) 
=indgen(512)*3.133    +1798 

t=100, 2000, 100  do  begin 
i-i+1 


Calculates  synthetic  Vegard-Kaplan 
spectra  at  temps  from  100  to  2000 
degrees  for  wavelengths  from  1800  to 
3400  A.   Results  convolved  with  MUSTANG 
instrument  slit  function 


int= 
intl 
int2 
int3 
int4 
int5 
int6 
int7 
int8 
int9 
intl 
intl 
intl 
intl 
intl 
intl 
intl 
intl 
intl 
intl 
int2 
int2 
int2 
int2 
int2 
int2 
int2 
int2 
int2 


wl*0 

=wl*0 

=wl*0 

=wl*0 

=wl*0 

=wl*0 

=wl*0 

=wl*0 

=wl*0 

=wl*0 

0=wl*0 

l=wl*0 

2=wl*0 

3=wl*0 

4=wl*0 

5=wl*0 

6=wl*0 

7=wl*0 

8=wl*0 

9=wl*0 

0=wl*0 

l=wl*0 

2=wl*0 

3=wl*0 

4=wl*0 

5=wl*0 

6=wl*0 

7=wl*0 

8=wl*0 


these  variables  provide  a  means  of  tweaking 
individual  transitions 


This  section  calculates  the  transitions 
into  the  total  spectrum. 

TRANSITION  WL 


band  systems,  summing  each  one 


vksyn, 
vksyn, 
vksyn, 
vksyn , 
vksyn, 
vksyn, 
vksyn, 
vksyn, 
vksyn, 
vksyn, 
vksyn, 
vksyn, 
vksyn, 
vksyn, 


t,6, 
t,4, 
t,13 
t,ll 
t,9, 
t,7, 
t,5, 
t,3, 
t,12 
t,10 
t,8, 
t,6, 
t,4, 
t,13 


l,wl, 

0,wl, 
,  5,wl 
,4,wl 
3,wl, 
2,wl, 
l,wl, 
0,wl, 
,  5,wl 
,4,wl 
3,wl, 
2,wl, 
l,wl, 
,6,wl 


int 

int 

,  int 

,  int 

int 

int 

int 

int 

,  int 

,  int 

int 

int 

int 

,  int 


1798 
1808 
1823 
1824 
1827 
1833 
1841 
1853 
1860 
1863 
1868 
1876 
1887 
1899 
81 


vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn , t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
int=int+ 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
int=int+ 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
int=int+ 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn , t , 
int=int+ 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
int=int+ 
vksyn , t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
int=int+ 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
int=int+ 


ll,5,wl 
2,0,wl, 
9,4,wl, 
7,3,wl, 
5,2,wl, 
3,l,wl, 
12,6,wl 
10,5,wl 
8,4,wl, 
1,0, wl, 
intl*l. 
6 , 3 ,wl , 
4,2,wl, 
13,7,wl 
ll,6,wl 
9,5,wl, 
2,1, wl, 
7,4,wl, 
0,0, wl, 
int2*2. 
5 , 3 ,wl , 
12,7, wl 
3,2,wl, 
10,6,wl 
8, 5,wl, 
1,1, wl, 
int3*l. 
6,4,wl, 
4,3,wl, 
13,8,wl 
ll,7,wl 
9,6,wl, 
2 , 2 , wl , 
7, 5,wl, 
5,4,wl, 
0,1, wl, 
int4*2. 
12,8,wl 
10,7,wl 
3 , 3 ,wl , 
8,6,wl, 
6 , 5 ,wl , 
l,2,wl, 
int5*l. 
13,9, wl 
4,4,wl, 
ll,8,wl 
9,7,wl, 
2, 3,wl, 
7,6,wl, 
5, 5,wl , 
12,9,wl 
0,2,wl, 
int6*2. 
10,8,wl 
3,4,wl, 
8,7,wl, 
6 , 6 ,wl , 
1, 3,wl, 
int7*l. 


,  int 

int 

int 

int 

int 

int 

,  int 

,  int 

int 

intl 

5 

int 

int 

,  int 

,  int 

int 

int 

int 

int2 

6 

int 

,  int 

int 

,  int 

int 

int3 

5 

int 

int 

,  int 

,  int 

int 

int 

int 

int 

int4 

6 

,  int 

,  int 

int 

int 

int 

int5 

5 

,  int 

int 

,  int 

int 

int 

int 

int 

,  int 

int6 

6 

,  int 

int 

int 

int 

int7 

5 


1900 
1901 
1905 
1912 
1923 
1936 
1939 
1943 
1950 
1954 

1960 
1973 
1980 
1983 
1989 
1990 
1998 
2010 

2011 
2024 
2027 
2029 
2038 
2047 

2050 
2065 
2067 
2071 
2079 
2085 
2090 
2106 
2109 

2115 
2122 
2125 
2133 
2147 
2148 

2161 
2166 
2167 
2177 
2189 
2191 
2209 
2213 
2216 

2222 
2231 
2236 
2253 
2258 

82 


vksyn, t , 
vksyn , t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn , t , 
vksyn , t , 
vksyn, t , 
vksyn, t , 
int=int+ 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
int=int+ 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
int=int+ 
vksyn, t , 
vksyn , t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
int=int+ 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 
int=int+ 
vksyn, t , 
vksyn, t , 
vksyn, t , 
int=int+ 
vksyn, t , 
vksyn, t , 
int=int+ 
vksyn, t , 
vksyn, t , 
int=int+ 
vksyn, t , 
int=int+ 
vksyn, t , 
vksyn, t , 
vksyn, t , 
vksyn, t , 


13 

,10,* 

l,int 

11 

,9,wl 

,  int 

4, 

5,wl, 

int 

9, 

8,wl, 

int 

7, 

7,wl, 

int 

2, 

4,wl, 

int 

12 

,10,  w 

1 ,  int 

5, 

6,wl, 

int 

10 

,9,wl 

,  int 

o, 

3,wl, 

int8 

in 

t8*2. 

6 

3, 

5,wl, 

int 

8, 

8,wl, 

int 

6, 

7,wl, 

int 

13 

,11,* 

l,int 

1, 

4,wl, 

int9 

in 

t9*l. 

5 

11 

,10,* 

1 ,  int 

4, 

6,wl, 

int 

9, 

9,wl, 

int 

7, 

8,wl, 

int 

2, 

5,wl, 

int 

12 

,11, w 

1 ,  int 

5, 

7,wl, 

int 

10 

,10,* 

1 ,  int 

0, 

4,wl, 

intlO 

in 

tl0*2 

.6 

8, 

9,wl, 

int 

3, 

6,wl, 

int 

13 

,12,* 

l,int 

6, 

8,wl, 

int 

11 

,ll,w 

1 ,  int 

1, 

5,wl, 

intll 

in 

tll*l 

.5 

9, 

10, wl 

,  int 

4, 

7,wl, 

int 

7, 

9,wl, 

int 

2, 

6,wl, 

int 

12 

,12,* 

1 ,  int 

5, 

8,wl, 

int 

10 

,11,* 

1,  int 

8, 

10, wl 

,  int 

0, 

5,wl, 

intl2 

in 

tl2*2 

.6 

3, 

7,wl, 

int 

13 

,13,* 

1,  int 

6, 

9,wl, 

intl3 

in 

tl3*. 

83 

11 

,12,* 

1 ,  int 

1. 

6,wl, 

intl4 

in 

tl4*l 

.5 

9, 

11, wl 

,  int 

4, 

8,wl, 

intl5 

in 

tl5*. 

67 

7, 

10, wl 

,intl6 

in 

tl6*. 

6 

12 

,13,* 

1,  int 

2, 

7,wl, 

int 

10 

,12,* 

1 ,  int 

5, 

9,wl, 

int 

2262 
2270 
2275 
2282 
2299 
2302 
2320 
2321 
2331 
2334 

2347 
2347 
2368 
2372 
2379 

2382 

2394 
2397 
2417 
2425 
2435 
2443 
2449 
2463 

2469 
2474 
2490 
2494 
2504 
2511 

2522 

2524 
2547 
2561 
2561 
2577 
2578 
2602 
2605 

2613 
2620 
2632 

2637 
2657 

2659 
2668 

2689 

2698 
2711 
2719 
2724 
83 


vksyn , 
vksyn , 
vksyn , 
int=in 
vksyn , 
vksyn , 
int=in 
vksyn , 
vksyn , 
vksyn, 
int=in 
vksyn, 
vksyn, 
vksyn, 
int=in 
vksyn, 
vksyn, 
int=in 
vksyn, 
int=in 
vksyn, 
int=in 
vksyn , 
vksyn, 
int=in 
vksyn, 
vksyn , 
vksyn , 
vksyn , 
vksyn, 
int=in 
vksyn , 
vksyn, 
vksyn, 
vksyn, 
vksyn , 
vksyn , 
vksyn , 
vksyn , 
vksyn , 
vksyn , 
vksyn , 
int=in 
vksyn , 
vksyn , 
vksyn, 
vksyn , 
vksyn , 
int=in 
vksyn , 
vksyn , 
vksyn, 
vksyn, 
vksyn, 
vksyn , 
vksyn, 
vksyn , 
vksyn, 
vksyn, 
int=in 


t,8,l 
t,13, 
t,0,6 
t+int 
t,3,8 
t,ll, 
t+int 
t,6,l 
t,9,l 
t,l,7 
t+int 
t,4,9 
t,7,l 
t,12, 
t+int 
t,10, 
t,2,8 
t  +  int 
t,5,l 
t+int 
t,8,l 
t+int 
t,13, 
t,0,7 
t+int 
t,3,9 
t,ll, 
t,6,l 
t,9,l 
t,l,8 
t+int 
t,4,l 
t,12, 
t,7,l 
t,10, 
t,2,9 
t,5,l 
t,13, 
t,8,l 
t,H, 
t,3,l 
t,0,8 
t+int 
t,6,l 
t,9,l 
t,12, 
t,4,l 
t,l,9 
t+int 
t,7,l 
t,10, 
t,2,l 
t,5,l 
t,13, 
t,8,l 
t,ll, 
t,3,l 
t,6,l 
t,0,9 
t+int 


l,wl, 
14, wl 
,wl,i 
17*2. 
,wl,  i 
13, wl 
18*. 7 
0,wl, 
2,wl, 
,wl,  i 
19*1. 
,wl,i 
l,wl, 
14, wl 
20*. 6 
13,  wl 
,wl,  i 
21*. 6 
0,wl, 
22*. 5 
2,wl, 
23*. 5 
15,wl 
,wl,  i 
24*2. 
,wl,  i 
14, wl 
l,wl, 
3,wl, 
,wl,  i 
25*1. 
0,wl, 
15, wl 
2,wl, 
14, wl 
,wl,  i 
l,wl, 
16, wl 
3,wl, 
15,wl 
0,wl, 
,wl,  i 
26*2. 
2,wl, 
4,wl, 
16, wl 
l,wl, 
,wl,  i 
27*1. 
3,wl, 
15, wl 
0,wl, 
2,wl, 
17, wl 
4,wl, 
16, wl 
l,wl, 
3,wl, 
,wl,  i 
28*2. 


int 

,  int 

ntl7 

6 

nt 

,intl8 

int 

int 

ntl9 

5 

nt 

int 

,int20 

7 

,  int 

nt21 

7 

int22 

int23 

,  int 

nt24 

6 

nt 

,  int 

int 

int 

nt25 

5 

int 

,  int 

int 

,  int 

nt 

int 

,  int 

int 

,  int 

int 

nt26 

6 

int 

int 

,  int 

int 

nt27 

5 

int 

,  int 

int 

int 

,  int 

int 

,  int 

int 

int 

nt28 

6 


;  274  8 
;  2761 
;  2762 

;  2767 
;  2782 

;  2783 
;  2810 
;  2818 

;  2826 
;  2845 
;  2848 

;  2875 
;  2877 

;  2887 

;  2909 

;  2916 
;  2937 

2938 
2942 
2951 
2976 
2998 

3002 
3013 
3017 
3046 
3062 
3068 
3087 
3087 
3119 
3128 
3133 

3137 
3159 
3196 
3198 
3200 

3210 
3235 
3269 
3270 
3276 
3285 
3315 
3342 
3345 
3353 

84 


vksyn, t , 9 , 15 ,wl , int 
vksyn, t, 12,17, wl , int 

vk( i  )=convol( int, slit ) 


3364 
3399 


endf or 

close , 1 

close , 3 

return 

end 

@[ synspecjvksyn 
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;Walden,1991 

PRO  NODLTCALC,WL, INT, SLIT 

OPENU ,10,' NODLT . DAT ' 
DLT=ASSOC( 10 , FLTARR( 512 ) ) 


WL=INDGEN( 512)*3.133  +  1798 


This  procedure  calculates  synthetic 
nitric  oxide  delta  band  spectra 
at  temps  from  100  to  2000  deg  and 
convolves  results  with  MUSTANG  instru- 
ment slit  function. 


OPENR, 11, ' [WALDEN] SLIT. DAT' 
SL=ASSOC( 11,FLTARR( 11 ) ) 
SLIT=SL(1 ) 


1  =  0 
FOR 

T=100, 2000, 100  DO  BEGIN 

1  =  1  +  1 

;  Wave] 

.e 

INT=WL*0 

DELSYN,T,0,0,WL,INT 

•  1909. 

8 

DELSYN,T,0,1,WL,INT 

•  1980. 

5 

DELSYN,T,0,2,WL, INT       , 

•  2055. 

9 

DELSYN,T,0, 3,WL, INT      , 

•  2134. 

8 

DELSYN,T,0,4,WL, INT 

■  2219. 

6 

DELSYN,T,0, 5,WL, INT       , 

-  2310. 

1 

DELSYN,T,0,6,WL, INT       , 

•  2407. 

6 

DELSYN,T,0,7,WL, INT      , 

■  2509 

DELSYN,T,0,8,WL,INT      , 

•  2620 

DELSYN,T,0,9,WL, INT       , 

•  2739 

DELSYN,T,0,10,WL, INT      , 

■  2866 

DELSYN,T, 0,11 ,WL, INT      , 

•  3003 

DELSYN,T,0,12,WL,INT     , 

•  3152 

DELSYN,T,0,13,WL, INT     , 

•  3313 

DELSYN,T, 0,14, WL, INT      , 

•  3488 

DLT( I )=CONVOL( INT, SLIT) 


ENDFOR 

CLOSE, 10 
CLOSE, 11 

RETURN 

:nd 

'[SYNSPEC]DELSYN 
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Walden,1991 

This  program  derives  the  synthetic  spectrum  using  tweaking  factors 

first  determined  by  Mack(1991)  for  1800-2100A,  and  extended  the  full  range 

of  the  MUSTANG  data. 

RO  NOGAMCALC,WL, INT, SLIT 


OPENU ,10,' NOGAM . DAT ' 
GAM=ASSOC(10,FLTARR( 512) ) 
WL=INDGEN( 512)*3.133  +  1798 


This  procedure  calculates  synthetic 
nitric  oxide  gamma  band  spectra 
at  temps  from  100  to  2000  deg  and 
convolves  results  with  MUSTANG  instru- 
ment slit  function. 


OPENR,ll, ' [WALDEN] SLIT. DAT' 
SL=ASSOC( 11 , FLTARR( 11 ) ) 
SLIT=SL(1) 


1  =  0 
FOR 


system. 


T=100, 2000, 100  DO  BEGIN 

1  =  1  +  1 

;These  variables  allow 

;individual  transition' 

INT=WL*0 

INTl=WL*0 

INT2=WL*0 

INT3=WL*0 

INT4=WL*0 

INT5=WL*0 

INT6=WL*0 

INT7=WL*0 

INT8=WL*0 

INT9=WL*0 

INT10=WL*0. 

INT11=WL*0. 

INT12=WL*0. 

INTl3=WL*0. 

INTl4=WL*0. 

INT15=WL*0. 

INT16=WL*0 . 

INT17=WL*0. 

INT18=WL*0. 

INT19=WL*0. 

INT20=WL*0. 

INT21=WL*0. 

INT22=WL*0. 

INT23=WL*0. 

This  section  calculates 

Transition 


for  the  adjusting 
s  probability 


of 


the  spectrum  from  each  transition's  band 
WL 


GAMSYN,T,3,0,WL,INT 

t 

1956 

1 

GAMSYN,T, 3,1,WL,INT 

t 

2029 

8 

GAMSYN , T , 2 , 0 , WL , INT9 

i 

2046 

5 

INT=INT+INT9* .75 

GAMSYN, T, 3,2,WL, INT 

i 

2108 

5 

GAMSYN , T , 2 , 1 , WL , INT 

} 

2128 

4 

GAMSYN, T,l, 0,WL, INTl 

? 

2148 

1 

INT=INT+INTl*1.125 

GAMSYN,T, 3, 3,WL,INT 

i 

2192 

8 

GAMSYN , T , 2 , 2 , WL , INT2 

i 

2215 

4 

87 


INT=INT+INT2 
GAMS YN,T, 1,1 
INT=INT+INT3 
GAMSYN,T,0, 0 
INT=INT+INT4 
GAMSYN,T, 3,4 
GAMSYN,T,2,3 
INT=INT+INT5 
GAMSYN,T,1 ,2 
INT=INT+INT1 
GAMS YN,T, 0,1 
INT=INT+INT6 
GAMSYN,T, 3,5 
GAMSYN,T,2,4 
GAMSYN,T,1 , 3 
INT=INT+INT7 
GAMSYN,T,0,2 
INT=INT+INT8 
GAMSYN,T, 3,6 
GAMSYN,T,2,5 
GAMSYN,T,1 , 4 
INT=INT+INTl 
GAMSYN,T,0,3 
INT=INT+INT1 
GAMSYN,T, 3,7 
GAMSYN,T,2,6 
GAMSYN,T,1 , 5 
INT=INT+INT1 
GAMSYN,T, 3,8 
GAMSYN,T,0,4 
INT=INT+INT1 
GAMSYN,T,2,7 
GAMS YN,T, 1,6 
INT=INT+INTl 
GAMSYN,T, 3,9 
GAMSYN,T,0, 5 
INT=INT+INT1 
GAMSYN,T,2,8 
GAMS YN,T, 1,7 
INT=INT+INT1 
GAMSYN,T, 3,1 
GAMS YN,T, 0,6 
INT=INT+INT1 
GAMSYN,T,2,9 
GAMSYN,T,1 ,8 
INT=INT+INTl 
GAMSYN,T,3,1 
GAMS YN,T, 0,7 
INT=INT+INT2 
GAMS YN,T, 2,1 
GAMS YN,T, 1,9 
INT=INT+INT2 
GAMS YN,T, 3,1 
GAMS YN,T, 2,1 
GAMS YN,T, 0,8 
INT=INT+INT2 
GAMSYN,T, 3,1 
GAMSYN,T,1, 1 
INT=INT+INT2 


*1.26 

,WL, INT3 

;  2238 

.3 

*1.3 

,WL, INT4 

;  2261 

.8 

*1.2 

,WL,INT 

;  2283 

.0 

,WL,INT5 

;  2308 

.5 

*1.6 

,WL, INT10 

;  2  3  35 

0*1.8 

,WL,INT6 

;  2362 

.2 

*1.2 

,WL, INT 

;  2378 

.8 

,WL, INT 

;  24  09 

,WL,INT7 

;  24  38 

.6 

*1.6 

,WL,INT8 

;  2470 

.1 

*1.2 

,WL, INT 

;  2481 

.2 

,WL, INT 

;  2516 

.4 

,WL, INT17 

;  2550 

.3 

7*1.3 

,WL, INTll 

;  2586 

.3 

1*1.2 

,WL,INT 

;  2591 

,WL, INT 

;  2629 

.9 

,WL, INT12 

;  2670 

1 

2*1.3 

,WL, INT 

;  2709 

,WL, INT13 

;  2712 

0 

3*1.2 

,WL,INT 

;  2753 

6 

,WL, INT14 

;  2800 

0 

4*1.3 

,WL, INT 

;  2836 

,WL, INT15 

;  28  48 

2 

5*1.2 

,WL,INT 

;  2888 

2 

,WL, INT16 

;  29  4  0 

5 

6*1.3 

0,WL, INT 

;  2973 

,WL, INT18 

;  2997. 

6 

8*1.2 

,WL,INT 

;  304  4. 

3 

,WL, INT19 

;  3112. 

4 

9*1.3 

1,WL, INT 

;  3120. 

6 

,WL, INT20 

;  3170. 

7 

0*1.2 

0,WL, INT 

;  3201. 

1 

,WL, INT21 

;  3278. 

5 

1*1.3 

2,WL,INT 

;  3  303. 

0 

1,WL,INT 

;  3361 

,WL,INT22 

;  3  375. 

5 

2*1.2 

3,WL, INT 

;  3  4  56 

0,WL, INT23 

;  3  4  58. 

5 

3*1.3 

88 


GAM( I )=CONVOL( INT, SLIT) 

ENDFOR 

CLOSE, 10 
CLOSE, 11 

RETURN 

END 

@[SYNSPEC]GAMSYN 
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;Walden,1991 
PRO  NOEPSCALC 


OPENU , 1 , ' NOEPS . DAT ' 
EPS=ASSOC( 1 , FLTARR( 512 ) ) 


WL=INDGEN( 512 )*3.133  +  1798 


This  procedure  calculates  synthetic 
nitric  oxide  epsilon  band  spectra 
at  temps  from  100  to  2000  deg  and 
convolves  results  with  MUSTANG  instr- 
ment  slit  function,  where  EPS(l) 
corresponds  to  T=100,etc.   It  uses 
scale  factors  determined  by  Bosserman 
except  for  (0,1),  which  Bosserman 
used  .84  to  tweak 


OPENR,2, ' [WALDEN]SLIT.DAT' 
SL=ASSOC ( 2 , FLTARR (11)) 
SLIT=SL(1 ) 


1  =  0 

FOR 

T=100, 2000, 100  DO  BEGIN 

1  =  1  +  1 

;  Wavele 

INT=WL*0 

INTl=WL*0 

INT2=WL*0 

INT3=WL*0 

INT4=WL*0 

EPSSYN,T,1 ,0,WL, INT 

;  1799. 

6 

EPSSYN,T,2,2,WL, INT 

;  18  49 

5 

EPSSYN,T,1,1,WL, INT 

;  1863. 

2 

EPSSYN,T,0,0,WL,INT1 

;  1876. 

6 

INT=INT+INT1*.63 

EPSSYN,T,2, 3,WL, INT 

;  1918 

EPSSYN,T,1,2,WL, INT 

;  1932. 

3 

EPSSYN,T,0,1,WL, INT2 

;  194  5. 

0 

INT=INT+INT2*.90 

EPSSYN,T,1, 3,WL, INT 

;  1998. 

6 

EPSSYN,T,0,2,WL,INT3 

;  2  017. 

5 

INT=INT+INT3*1 .03 

EPSSYN,T,1,4,WL,INT 

;  2  073. 

0 

EPSSYN,T,0, 3,WL, INT4 

;  2094. 

5 

INT=INT+INT4*1 .11 

EPSSYN,T,1, 5,WL,INT 

;  2157. 

5 

EPSSYN,T,0,4,WL,INT 

;  2176. 

1 

EPSSYN,T,2,8,WL, INT 

;  2294. 

0 

EPSSYN,T,2,9,WL,INT 

;  2  385. 

0 

EPS( I )=CONVOL( INT, SLIT) 


ENDFOR 

CLOSE, 2 
CLOSE, 1 

ETURN 

ND 

[ synspec]EPSSYN 
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;Walden,1991 

pro  01  2972calc  ;This  program  creates  a  synthetic  model  for  a 

;  neutral  atomic  oxygen  emission 

openu,l, 'OI_297  2.dat' 
dat=assoc(l,fltarr(512)  ) 

openr ,2, ' [waldenjslit.dat' 

sl=assoc(2,fltarr(ll)  ) 

slit=sl(4)  ;  adjusts  wl  to  2972  vice  2972.875 

hold=fltarr(512) 

hold=indgen( 512)* 0.0 
hold(375)=1.0 

dat( 0 )=convol ( hold, slit )/3 . 133   ; spectra=l/wl  bin  so  that  scale  factor  is 

;absolute  intensity 

close ,  2 
close ,  1 

return 
end 
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;Walden,1991 
pro  OII_2470calc 

openu,l , 'OII_2470.dat' 
dat=assoc(l,fltarr(512)  ) 

openr , 2 , ' [ walden ] sli t .dat ' 
sl=assoc(2,fltarr(ll) ) 
slit-sl(l) 


;program  calculates  synthetic  spectrum 
;for  0+  atomic  emission  @  2470A 


;this  makes  wl  centered  @  2968. 5A 


hold=fltarr(512) 

hold=indgen( 512)*0.0 
hold(214)-1.0 

dat (0)=convol( hold, slit)/3. 133 


;spectra  is  1/wl  bin  so  scale  factor 
;is  absolute  intensity 


close , 2 
close , 1 

return 
end 
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;Walden,1991 

;This  program  creates  a  syntheic  model  for  the  0++  atomic  emission  @  2853A. 

pro  OIII_2853calc 

openu,l, 'OIII_2853.dat' 
dat=assoc(l,fltarr(512) ) 

openr ,  2 , ' [waldenjslit.dat' 
sl=assoc(2,fltarr(ll) ) 
slit-sl(l) 

hold=fltarr(512) 

hold=indgen( 512)*0.0 
hold(337)=1.0 

spike=convol( hold, slit) 

dat( 0 )=spike/total ( spike )/3 . 133  ;dat(0)  is  regular  convolved  spike 

spike ( 338 : 510 )=spike ( 339 : 511 )    ;dat(l)  is  spike  with  adjustment  for  apparent 

;  dropped  point  in  data. 
dat(l)=( spike/total (spike) )/3.133 

close ,  2 
close ,  1 

return 

end 
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;Walden,1991 

;This  program  calculates  a  synthetic  spectrum  for  the  Nil  doublet 
;with  the  ratio  based  on  the  findings  of  Bucsela  and  Sharp  (1989) 
pro  NII_doubcalc 

openu,l, 'NII_2143doub.dat' 

nii=assoc(l,fltarr(512)  ) 

openr ,  2 , ' [walden]slit.daf 
sl=assoc(2,fltarr(ll)  ) 

hold=fltarr(512) 

hold(109)=.58  ;2139.68A    line 

intl=convol ( hold, si ( 1 ) )  ;wl  is  in  center  of  pixel 

hold2=fltarr(512) 

hold2(110)-1.0  ;2143.55  line 

int2  =  convol  (  hold2  ,  si  (  3  )  )  ;wl  is  1/4  to  the  ri.ght  of  center  of  pixel 

int=intl+int2 

nii(l)-(int/total(int) )/3.133 

;total  of  spectral  area  is  1/wl  bin,  so 
;scale  factor  will  be  the  absolute  intensity 

close , 2 
close , 1 
return 
end 
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;Walden,1991 

;This  program  creates  a  synthetic  spectrum  of  atomic  emissions  observed  in 

;the  March  1990  MUSTANG  data. 

pro  otherspec 

openr , 1, ' [waldenjslit.dat' 
slit=assoc(l,fltarr(ll) ) 

openu, 2, 'otherspec. dat' 
spect=assoc(2,fltarr(512) ) 

spec  =  f ltarr ( 512  ) 

numlines=20  ;number  of  contributing  atomics 

;wavelengths  determined  from  Bureau  of  Standards 
;listing  of  atomic  emissions  and  NASA  listing  of  atomic 
; emissions 


Wavelength 

Species 

Pixel 

1846.0 

NIII 

15 

2182.6 

Oil 

123 

2190.6 

Oil 

125 

2337.8 

ARII 

172 

2384.6 

OIV 

187 

2445.6 

Oil 

207 

2490.3 

Nil 

221 

2509.9 

Nil 

227 

2665.8 

OIII 

277 

2709.8 

Nil 

291 

2741.3 

Oil 

301 

2885.3 

Nil 

347 

2973.6, 

2974 

.6, 

Nil  multiplet 

375-377 

2975.7, 

2975 

.9, 

2976.9, 

2977 

.3, 

2978.4, 

2978 

.6 

pix-[ 15 ,123, 125, 172, 187, 207, 221, 227, 277 ,291, 301, 347,$ 

37  5,375,376,376,376,376,377,377] 

tht-[ 10, 2, 1,10, 9, 3, 3, 4, 5, 6, 4, 4, 2, 3, 3, 4, 4, 4, 5, 5] 
3lno=[ 3, 4, 3, 3, 4, 4, 1,3, 1,1, 1,1, 3, 2, 4, 1,3, 2, 4, 4] 

for    i=0 , numlines-1    do   begin 

blank=fltarr(512) 

blank(pix(i) )=.l*int(i) 

impulse=convol (blank, slit(slno(i) ) ) 

spec=spec+ impulse 
?ndf or 

pect(0)=indgen( 512 )*3. 13 3+17  98 
;pect( 1 )=( spec/total ( spec ) )/3 . 133        ;so  that  scale  factor  is  absolute  int 

:lose ,  2 
lose ,  1 
eturn 
nd 
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APPENDIX  C 

Appendix  C  is  the  programs  used  to  determine  the  standard 
deviations  of  the  means (sdoms)  for  the  data.  The  sdoms  are 
used  to  determine  the  weighting  factors  for  the  fitting 
routine  by  the  program  titled  Mweights2,"  which  is  also 
included. 
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;Walden,1991 

pro  m90sdev ,alt , indexl , index2 

;This  procedure  calculates  the  std.dev  and  std  dev  of  the  mean 

;for  the  averaged  spectra. 

;Inputs  consist  of  the  requested  altitude  (1-25)  and  its  corresponding 

;beginning  and  ending  index  numbers. 

;Outputs  are  the  SDOM  (unshifted  and  shifted)  and  are  written  to 

;um90_sdom.dat  and  m90_sdom.dat,  respectively 

iavg=0 
5t='  ' 

sigma=dblar r ( 512 ) 
5igmam=dblarr ( 512 ) 

Dpenr , 1 , 'dubO : [ clay ton  Jmustang_raw.dat ' 

3at=assoc( 1 , intar r ( 512 ) )  ;opens  the  raw  data  file 

Dpenr , 2 , ' [ mustang]um90_avg .dat '  ;opens  the  averaged  spectrum  file 

avg=assoc( 2 , f ltarr ( 512 ) )  ;  (unshifted  version) 

Dpenu, 3 ,' [mustang Jum90_sdm.dat '  ;opens  the  file  to  write  SDOM  into 

sdom=assoc( 3,fltarr(512)) 

i=index2-indexl 
Eor  j  =  0,n  do  begin 

print , indexl+ j 

print, 'Do  you  want  to  include  this  record' 

print, 'Type  <CR>  for  yes,  N  for  no.' 

read, st 

if  st  eq  'n'  then  goto,labell 

sigma=sigma  +  ( dat ( indexl+ j )-avg( alt ) ) "2 

navg=navg+l 

labell :dummy=0 
?ndf  or 
:ree=navg-l 

;igma=sqrt ( sigma/f ree ) 
igmam=sigma/sqrt ( navg ) 
dom( alt )=sigmam 

lose, 1 
lose , 2 
lose ,  3 

eturn 
nd 
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;Walden,1991 
pro  weights2 

;This  program  checks  for  values  where  the  detector  saturated,  and  gives 
;weights  based  on  the  deviations  of  the  mean  being  equal  to  the  data  at 
;those  points.   It  also  gives  a  small  weight  factor  to  those  elements 
/corresponding  to  nitrogen  second  positive  wavelengths  (sets  sdom  equal  to 
/largest  data  value  in  N2  2nd  pos  bands).  Since  sdom  was  given  value  of  data  at 
/dropped  points,  it  gives  small  weights  to  dropped  values  (with  the  exception 
;of  2972A) 

openr ,1, ' [mustang] m9  0_intsdm.dat' 
dom=assoc( l,fltarr(512)) 

openr, 2, ' [mustangjm90_int.dat' 
dat=assoc(2,fltarr(512)  ) 

openu, 3 , 'm9  0_domwt . dat ' 
domwt=assoc( 3,fltarr(512)) 


openr ,  4 , ' [mack ]mack_cal 
ss=assoc(4,fltarr(512)  ) 
sens=ss( 22 ) 


dat 


openr , 5 , ' [mustang ] m9  0_avg.dat ' 
dd=assoc(5,fltarr( 512)  ) 

big=1800 

holddom=fltarr( 512) 
holddomwt=f ltarr( 512) 
holddat=fltarr(512) 
sat=900.00/sens 


;approx  N2  2nd  pos  largest  datapoint  value 


for  i 

=1,25  do  begin 
holddom=dom( i ) 
holddat=dat( i ) 

for  j=0,511  do 

begin 

;if  detector 

was  saturated, 

make 

sdom  equal  to 

data  value 

if 

(holddat 

j)  gt  sat(j)" 

I  then 

holddom( - 

j)-holddat(j) 

;if  wavelength  is  ] 

.n  N2  2nd 

pos 

r 

make  sdom  equal  to  big 

wavel 

ength 

if  (] 

eq 

277 

1  or  (  : 

eq 

278) 

or 

eq 

279 

)  then 

holddom( ; 

)=big 

r  2669 

if  (: 

eq 

318 

1  or  (  ; 

eq 

319) 

or 

eq 

320 

)  then 

holddom( ; 

)  =big 

;2796 

if  (: 

eq 

322 

l  or  (  ■ 

1  eq 

323) 

or 

eq 

324 

)  then 

holddom( ; 

)  =big 

;2811 

if  (: 

eq 

362 

'  or  (  - 

eq 

363) 

or 

eq 

364 

I  then 

holddom( ■ 

)=big 

;2934 

If  (; 

eq 

368 

I  or  (  ; 

eq 

369) 

or 

eq 

370 

then 

holddom( ; 

)=big 

•2953 

Lf  (] 

eq 

374 

or  (  - 

1  eq 

375) 

or 

eq 

376 

\    then 

holddom( ; 

)=big    i 

•2973 

if  (: 

eq 

409 

or  (  ; 

1  eq 

410) 

or 

eq 

411 

then 

holddom( ; 

)=big    j 

•3083 

if  (: 

eq 

417 

1  or  (  ; 

1  eq 

418) 

or 

eq 

419 

then 

holddom( ; 

)=big    , 

•3107 

-f  (: 

eq 

425 

1  or  (; 

1  eq 

426) 

or 

eq 

427 

then 

holddom( ; 

)=big    j 

•3132 

-f  (- 

eq 

433 

'  or  ( - 

1  eq 

434) 

or 

eq 

435 

then 

holddom( - 

)=big    ] 

•3157 

f  (- 

eq 

461 

1  or  (; 

1  eq 

462) 

or 

eq 

463 

then 

holddom( \ 

)=big    ; 

3244 

f  (- 

eq 

470 

1  or  (; 

eq 

471) 

or 

eq 

472 

then 

holddom( ; 

)=big    ( 

3274 

f  (; 

eq 

480 

l  or  ( - 

eq 

481) 

or 

- 

eq 

482 

then 

holddom( i 

)=big    ; 

3305 

f  (: 

eq 

490 

•  or  ( ; 

eq 

491) 

or 

: 

eq 

492 

then 

holddom( ; 

)=big    ; 

3337 

f  (- 

eq 

501 

or  ( ; 

eq 

502) 

or 

: 

eq 

503, 

then 

holddom( ; 

)=big    j 

3370 

f  (j  gt  493)  then  holddom( j )=big 

now  give  weights  their  values 

if  holddom(j)  eq  0.  then  holddomwt ( j ) =1 

98 


;last  spike 


if  holddom(j)  ne  0.  then  holddomwt ( j ) =1 ./holddom( j ) "2 
endf or 
Correct  for  where  low  value  of  data  point  at  dropped  word  gives  false 
high  weight  to  point  by  giving  dropped  points  1/10  the  maximum  weight 
of  any  other  point 

(except  for  2972A  -pixels  371  to  376  and  2850A  -333  to  336). 
dropped=where ( dd( i )  eq  -1.) 
ho ld= holddomwt 
hold(dropped)=0.0 
holddomwt ( dropped )=.l*max(hold) 
holddomwt (371: 376) -max (hold) 
holddomwt ( 3  33 : 3  36 )=max ( hold ) 
assign  weight  to  file  m90_domwt.dat 
domwt ( i )=holddomwt 
endf  or 
close ,  1 
close ,  2 
close ,  3 
close ,  4 
close ,  5 

return 
end 
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;Walden,1991 

;This  program  shifts  the  MUSTANG  data  array  according  to  dropped  pixels,  as 

;input  in  the  procedure  shiftarray. 

pro  shiftdata 


openr , 1 , 
mm=assoc 

openu , 2 , 
dd=assoc 

openr , 3 , 
uu=assoc 

openu, 4 , 
ss=assoc 


[mustang ]um9  0_avg . dat ' 
l,fltarr(512)  ) 

[mustang ] m9  0_avg.dat ' 
2,fltarr(512)  ) 

[ mustang  Jum9  0_sdm.dat ' 
3,fltarr(512)  ) 

[mustang ] m9  0_sdm.dat ' 
4,fltarr(512)  ) 


wl=indgen( 512 )*3. 133+1798 

dd(0)=wl 
ss( 0)=wl 

for  i=l,25  do  begin 

shiftarray,mm(i) ,holdl 

dd(i)=holdl 

shiftarray, uu(i) ,hold2 

ss( i )=hold2 
endf or 

close ,  4 
close ,  3 
close ,  2 
close ,  1 

return 
and 
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Walden,1991 

This  program  shifts  an  array  of  data  according  to  the  "drop"  array  of 
pixels . 
pro  shif tar  ray , inar r , outar r 

drop=[ 462, 4 46 ,430, 414 ,398, 382 ,366, 3 50 ,3 50, 3 34, 318, 3 02, 286, 27 0,254, 2 39, 22 3, 223,$ 
207,191,175,159,143,127,111,95,95,80,80,64,48,32,16] 


ndrops=n_elements ( drop ) 
outarr=inar r 
length=n_elements( inar r  ) 

for  i=0,ndrops-l  do  begin 

;        This  section  is  to  put  a  -1  in  dropped  values  and  shift  as  required 

outarr(drop(i)+l)=outarr(drop(i):length-2) 

outarr(drop( i ) )=-1.0 
endf or 

return 
end 
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;Walden,1991 

;This  program  takes  the  MUSTANG  data  from  the  averaged,  shifted  records 
;and  puts  in  the  appropriate  average,  and  adjusts  for  instrument  response 
pro  average 

openr,l, ' [mustang ]m9  0_avg . dat ' 
mm=assoc (l,fltarr(512)  ) 

openr,2, ' [mustang ] m9  0_sdm.dat ' 
ss=assoc(2,fltarr(512)  ) 

openu, 3 , ' [mustang ]m9  0_int .dat' 
dat=assoc(3,fltarr(512)  ) 

openu, 4, ' [ mustang ]m90_intsdm. dat ' 
dev=assoc( 4 , f ltarr ( 512  )  ) 

openr ,  5 , ' [mack ]mack_cal .dat ' 

sen=assoc(5,fltarr(512)  ) 

sens=sen( 22  ) 
boolean=f ltarr ( 512  ) 
dat(0)=mm(0) 
dev(0)=ss( 0) 

for  i=l,25  do  begin 

;get  array  of  where  dropped  words  are 

points=mm( i ) 

index=where ( points  eq  -1.0) 
jcorrect  for  sensitivity 

points=points/sens 
;assign  average  for  dropped  data  values 

points ( index )=( points ( index+1 )+points( index-1 ) )/2 . 
;correct  averages  for  double  drops 

boolean( index)=l 

double=where (boolean*shift(boolean,-l) ) 

points ( double )=(2*points( double-1 )+points ( double+2 ) )/3 . 

points  (  double+1 )  =  ( points ( double-1 )+2*points( double+2 ) )/3 . 
:assign  value  of  data  to  dropped  word  sdom's 

errors=ss ( i ) 

errors! index)=points( index) 
zero  first  six  elements 

points( 0:5)=0.0 

errorst 0: 5)«0.0 
assign  data  to  m90_int.dat 

dat ( i )=points 
correct  for  sensitivity  and  assign  sdom  to  m90_intsdm.dat 

dev( i )=errors/sens 
ndf  or 

:lose ,  5 
lose ,  4 
lose ,  3 
lose ,  2 
lose ,  1 

eturn 
nd 
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;Walden,1991 
pro  units 

openr , 1 , ' [mustanglm90_avg.dat' 
a=assoc(l,fltarr(512)  ) 

openr , 4 , ' [ mustang jm9  0_sdm.dat ' 
c=assoc(4,fltarr(512)  ) 

openr , 2 ,'[ mack ]mack  cal.dat' 
b=assoc(2,fltarr(512)  ) 
sens=b( 22 ) 


;reads  out  avg  spectra,  shifted 

; reads  out  SDOMs,  shifted 
;calibration  to  correct  units  to  intensity 


openw, 3, ' [mustangjm9  0_int.dat' ,512*4 
int=assoc(3,fltarr(512)  ) 

openw, 5 , ' [mustang ]m9  0_intsdm. dat ', 512*4 
intdom=assoc( 5,fltarr(512)) 

for  i=l,25  do  begin 

int ( i )  =  a ( i )/sens 

intdom(i)=c(i)/sens 
endf or 

close ,  1 
close ,  2 
close ,  3 
close ,  4 
rlose , 5 

return 
snd 


;file    to    accept    corrected   values 
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;Walden,1991 

;This  program  assigns  zeros  to  the  first  six  values  of  an  array  of  512 

pro  zerosix , f name 

openu, 1 , f name 
a=assoc(l,fltarr(512) ) 

temp=fltarr( 512) 

for  i=l,25  do  begin 

temp=a ( i ) 

temp( 0:5)=0.0 

a( i )=temp 
endf or 

close ,  1 

return 

end 
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APPENDIX  D 

Appendix  D  shows  the  programs  used  to  convert  the  scale 
factors  obtained  from  the  fitting  routine  to  absolute 
intensities.   The  eguations  were  derived  from  Barth  (1965). 

The  equation  used  for  converting  the  scale  factor  to 
intensity  for  the  VK(0,5)  band  is: 

I05-VKCALC(0,5)  *SF{alt)  , 

where  I0 5  is  the  band  intensity  in  Rayleighs,  VKCALC(0,5)  is 
the  total  intensity  of  the  synthetic  VK(0,5)  band  in 
Rayleighs,  and  SF(alt)  is  the  VK  scale  factor  for  the  altitude 
desired.  The  factor  VKCALC  is  just  the  area  under  the 
synthetic  spectrum  curve  for  the  transition  and  is  obtained 
using  the  IDL  command  TOTAL ()  and  multiplying  by  wavelength 
bin  size,  as  seen  in  the  program  VK_ABSINT. 

The  conversion  for  the  LBH  scale  factor  included  all 
bands,  and  is  performed  in  the  program  LBH_ABSINT  and  the 
subroutine  INTENCALC.   The  basic  equation  for  this  program, 
derived  from  equations  (4)  &  (8)  of  Barth  (1965)  is  given  by: 


LBHCALC(j,k)    Y,    Qj.v"vj.v* 
allv" 


■^■LBH"  .     3 


vjik-qjt0qjtk 

where  ILBH  is  the  total  intensity  of  LBH  emissions  in 
Rayleighs,  LBHCALC(j,k)  is  the  area  under  the  synthetic 
spectrum  curve  for  a  transition  band  (j,k)  ,  q.  .  is  the  Franck- 

J  /  K 
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Condon  factor  for  the  (j,k)  band,  and  v  is  the  corresponding 
frequency  for  each  transition.  The  summation  over  v"  is  to 
include  all  possible  upper  states.  The  result  of  this 
equation  was  essentially  the  same  for  any  (j,k).  This  thesis 
utilized  the  (4,12)  transition. 
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;Walden,1991 

;This  program  converts  the 

;an  absolute  intensity  for 

;Stored  in  vk_absint.dat,  the  values 

;the  (0,6)  in  record  3,  and  the  (1,8 

pro  vk_absint 

openr , 10 ,' final para .dat ' 

p=assoc(10,fltarr(25)  ) 

temp=p( 0 ) 

alt=p(21) 

vksf=p(13) 

delvksf=p(14) 

openu, 11 , ' vk_absint . dat ' 
dat=assoc(ll,fltarr(25)  ) 
dat(0)=alt 

wl=indgen( 512 )*3. 13 3+1798 

hold=fltarr(25) 

errhold=fltarr(25) 


scale  factor  from  the  nimpfit_all  routine  into 
a  particular  transition  band  of  Vegard-Kaplan. 

are  for  the  (0,5)  transition  in  record  1, 
in  record  5. 


for  ru 


n=2 , 6 
if 
vl  = 
v2  = 
ql  = 
q2  = 
twk 
end 
if 
vl  = 
v2  = 
ql  = 
q2  = 
twk 
end 
if 
vl  = 
v2  = 
ql  = 
q2  = 
twk 
end 
for 


endf or 


lose , 10 
lose ,  11 


+2  do  begin 

run  eq  2)  then  begin 


.9e-4 

.69e-l 

2.6 

f 

run  eq  4 )  then  begin 


;F-C  factors  are  not  currently  used 
;in  the  determination  of  intensity, 
;but  are  included  for  future  use 

F-C  factor  for  (0,0) 

F-C  factor  for  (0,5) 


.9e-4 

.89e-l 

2.6 

f 

run  eq  6) 


;F-C  factor 
;F-C  factor 


then  begin 


.32e 

.87e 

1.5 

f 

a=0, 


for  (0,0) 
for  (0,6) 


;F-C  factor  for  (1,0) 
;F-C  factor  for  (1,8) 


end 
dat 
dat 


or 

run- 
run) 


24  do  begin 

int=indgen( 512 )*0.0 

vksyn,temp(a) ,vl,v2,wl,int 

int=int*twk 

hold(a)=total( int ) *3 . 133*vksf ( a ) 

errhold(a)=total( int ) *3 . 133*delvksf ( a ) 

l)=hold 
=errhold 


eturn 

nd 

[ synspec ] vksyn 
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;Walden,1991 
pro  lbh_absint 
;  This  program 


determines  the  absolute  intensity  of  LBH  emissions  based 


on  the  fits  of  intensity  to  three  different  transition  bands  (the 


/basically  equal  in  Walden's  analysis) 
;is  based  on  th  (2,9),  record  3  is  the 

openr ,10, 'final para. da t ' 

p=assoc(10,fltarr(25)  ) 

temp=p( 0  ) 

lbhsf=p(9) 

dellbhsf=p(10) 

alt=p(21) 


Stored  in  record  1  of  lbh 
(4,12)  ,  record  5  is  ( 3,10)  ." 


results  were 
absint .dat 


openr , 11 , ' lbh_tots .dat' 
t=assoc(ll,fltarr(20) ) 
lbhtotint=t( 0) 

openu, 12 , ' lbh_absint .dat ' 
i=assoc(12,fltarr(25) ) 
i(0)=alt 

hold=fltarr(25) 
errhold=fltarr(25) 

for  run=l,3  do  begin 

if  (run  eq  1)  then  begin 

vl  =  2 

v2  =  9 

endif 

if  (run  eq  2)  then  begin 

vl  =  4 

v2  =  12 

endif 

if  (run  eq  3)  then  begin 

vl  =  3 

v2  =  10 

endif 

for  a=0,24  do  begin 

intencalc,temp(a) ,vl,v2,bandint 

intlbh=bandint*3.133*lbhsf ( a )/lbhtotint ( ( temp( a )/10  0 )-l ) 

errintlbh=bandint*3.133*dellbhsf ( a )/lbhtotint ( ( temp( a )/100 ) -1 ) 

hold(a)=intlbh 

errhold(a)=errintlbh 
endf or 


ndf  or 


i (2*run-l)=hold 
i(2*run)=errhold 


;put  absolute  intensities  in  appropriate  record 
;put  errors  in  correct  record 


lose , 10 
1  o  s  e  ,  1 1 
lose ,  12 

eturn 

nd 

[ synspec ] lbhsyn 
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;Walden,1991 

;This  procedure  is  called  by  the  lbh_absint  program  in  order  to  determine 
;the  intensity  under  a  particular  "hump"  (a  transition  band)  in  the  LBH  curve 
;The  inputs  vl  and  v2  are  the  upper  and  lower  states'  vibrational  levels 
;for  the  transition  to  be  determined.   Intlbh  returns  the  intensity  value, 
pro  intencalc , T,vl ,v2 , intlbh 

openr , 1 , ' [ synspec ] frnkcondn.dat ' 
fc=assoc(l,fltarr(15,30)  ) 
lbhfc=fc(0) 

openr, 2, ' [walden] lbhtrans_wl . dat ' 
wl=assoc(2,f ltarr( 15,  30)  ) 
lbhwl=wl(0) 

summa=0 . 0 

wl=indgen(512)*3.13  3+1798 
int=fltarr(512) 
lbhfreq=fltarr(15,30) 

for  3=0,14  do  begin 

for  k=0,29  do  begin 

if  (lbhwl(j,k)  ne  0)  then  lbhf req( j , k )=2 . 997e8/lbhwl ( j , k ) 

if  (lbhwl(j,k)  eq  0)  then  lbhf req( j , k ) =0  .  0 

endf or 
endf or 

for  i=0,29  do  begin 

summa=summa+ ( lbhfreq(vl,i ) "3)*lbhfc(vl,i ) 
endf or 

lbhsyn,T,vl,v2,wl,int 

intvlv2=total(int) 

intlbh=intvlv2*summa/(lbhf req(vl,v2) ~3*lbhfc(vl,0)*lbhfc(vl,v2)) 

close ,  1 
close ,  2 
i  return 
end 
@ [ synspec ] lbhsyn 
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APPENDIX  E 

Appendix  E  contains  plots  of  the  final  fits  obtained, 
along  with  the  program  used  to  print  them.  The  plots  are 
divided  into  40oA  increments,  from  I8O0A  to  3400A,  and  are  in 
order  of  increasing  altitude.  Each  plot  shows  the  fit  versus 
data  in  the  upper  view,  and  a  breakdown  of  the  fit  by  major 
contributors  in  the  second  view. 
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;Walden,1991 

pro  plotapp , al t , xmin , xmax 


openr , 1 , 'finalpara.dat' 
para=assoc(l,fltarr(25) ) 

height=para( 21 ) 

openr ,  2 , ' [mustang] m9  0_intnh.dat' 
d=assoc(2,fltarr(512) ) 

openr , 3 , ' f inal twk . dat ' 
f=assoc(3,fltarr( 512) ) 

temp=para( 0 ) 

it=( temp(alt-l) )/100 

openr , 4 , 'allvk.dat' 
vk=assoc(4,fltarr(512) ) 
pvk=para( 13 ) 

openr , 5 , 'all_lbh.dat' 
lbh=assoc (5,fltarr(512) ) 
plbh=para( 9) 

openr , 6 , 'nod It. dat' 
dlt=assoc( 6,fltarr(512)) 

openr, 7, 'noeps.dat' 
eps=assoc(7,fltarr(512) ) 

openr , 8 , ' no gamma ck .dat' 
gam=assoc(8,fltarr(512) ) 

pno=para( 11 ) 

openr ,9, 'OIII_2853.dat' 
oiii=assoc(9,fltarr(512) ) 
poiii=para ( 7 ) 

openr ,10, 'OII_2470.dat' 
oii=assoc(10,fltarr(512) ) 
poii=para ( 5 ) 

openr, 11, '01  2972.dat' 
oi=assoc(ll,Tltarr(512) ) 
poi=para( 3 ) 

openr, 12, 'NII_214  3doub.dat' 
n2=assoc(12,fltarr(512) ) 
pn2=para ( 15 ) 

wl=indgen(512)*3.13  3+179  8.00  0 


program  plots  data  vs.  fit  and  fit  by 
parts  in  format  used  by  Walden  thesis 
for  appendix,  alt  can  be  1  to  25,  xmin 
and  xmax  are  presumed  to  be  1800-2200, 
2200-2600,2600-3000,  or  3000-3400,  so 
other  combinations  would  need  adjust- 
ments to  the  plot  key 


;data  is  in  records  d=l  to  25 


;fits  are  in  records  f=l  to  25 


111 


set  plot , ' hp' 

device , /portrai t ,/plotter_on_of f 

device, xoffset=l. 27 

device ,yoffset=.635 

device ,xsize=17 .78 

device,ysize=24.13 

;TOP  OF  PAGE — DATA  AND  FIT 
set_xy , xmin , xmax 
set_viewport, .18, .82, .59, .91 

if  (alt  ge  18)  then  $ 

!mtitle='Data  &  Fit  at  ' +st rt r im(  ( f ix ( height ( al t-1 ) ) ) , 2 )  +  ' km  (downleg)' 

if  (alt  It  18)  then  $ 

lmtitle='Data  &  Fit  at  ' +strtrim( ( f ix( height ( alt-1 ) ) ) , 2 )+ ' km  (upleg)' 

!xtitle='Wavelength(A) ' 
iy title-' Intensity(R/A) ' 

dat=d(alt) 
fit=f (alt) 

1  linetype  =  0 
plot , wl ,dat 

!  linetype  =  2 
oplot , wl , fit 

;KEY  FOR  TOP  PLOT--xrange  must  be  determined  in  terms  of  pixels  so  that 
;yrange  can  be  known  for  placing  the  key  in  the  correct  place  on  the  graph 
;case  section  does  this  for  the  ranges  designated 

case  1  of 

(xmin  eq  1800)  and  (xmax  eq  2200):  begin 

ldat=dat(0:129) 

lfit=fit(0:129) 

end 
(xmin  eq  2200)  and  (xmax  eq  2600):  begin 

ldat=dat(129:257) 

lfit-fit(129:257) 

end 
(xmin  eq  2600)  and  (xmax  eq  3000):  begin 

ldat=dat(257:384) 

lfit=fit(257:384) 

end 
(xmin  eq  3000)  and  (xmax  eq  3400):  begin 

ldat=dat(384:511) 

lfit=fit(384:511) 

end 
else:  begin  ;key  will  probably  not  be  correct 

ldat=dat  ;this  section  must  be  manually  altered 

lfit=fit  ;to  concur  with  the  wl  range 

end 
endcase 

;this  section  designates  x  and  y  coordinates  for  the  sample  lines  in  the  key 

xcoordl=[xmin+20,xmin+50] 

xcoord2=[xmin+20,xmin+100] 

xcoord3=[xmin+120,xmin+200] 

ycoordl=[ . 96*max ( ldat ) , .96*max(ldat) ] 
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ycoord2= 
ycoord3= 
ycoord4= 
ycoord5= 
ycoord6= 
ycoord7= 


90*max( ldat ) 
84*max( ldat) 
78*max( ldat) 
72*max( ldat) 
66*max ( ldat ) 
60*max( ldat) 


.90*max( ldat) 
.84*max( ldat) 
.78*max(ldat ) 
.72*max( ldat) 
. 66*max ( ldat ) 
.60*max( ldat) 


! 1 ine type=0 

oplot,xcoordl,ycoordl 

! linetype=2 

oplot,xcoordl,ycoord2 

! Iinetype=0 

xyouts ,xmin+60 , . 995*ycoordl ( 0 ) 

xyouts ,xmin+60 , . 995*ycoord2 ( 0 ) 

;BOTTOM  OF  PAGE — PLOT  OF  PARTS  AND  FIT 


'Data' 
'Fit' 


set_xy , xmin,xmax 
set_viewport, .18, .82, .09,  .41 

!mti tle=' Data  and  Fit  by  Parts' 
lxtitle='Wavelength(A) ' 
!ytitle=' Intensi ty ( R/A) ' 

1  linetype  =  0 

plot ,wl ,dat 

oplot ,wl ,dat 

I linetype=l 

oplot, wl,vk( it) *pvk(alt-l) 

!  linetype  =  2 

oplot, wlflbh(it)*plbh(alt-l) 

!  linetype  =  3 

nitric=pno(alt-l )*(gam( it)+eps(it)+.25*dlt(it))/3.133e6 

oplot ,wl , ni trie  ;dash-dot 

!  linetype  =  4 

oplot, wl,oiii(0)*poiii(alt-l)  ;  solid 

I  linetype  =  5 

oplot, wl,n2(l)*pn2(alt-l)  ;  solid 

!  linetype  =  5 

oplot, wl,oii(0)*poii(alt-l)  ;  solid 

!  linetype  =  5 

oplot, wl,oi(0)*poi(alt-l)  ;  solid 

;key  for  bottom  plot--x  and  y  ranges  from  top  plot  used  for  locating  key 


;makes  data  a  bold  line 

;  d  o  1 1  e  d 

;dashed 


case 
(  xmi 


1  of 
n  eq 


1800)  and  ( xmax  eq  2200):  begin 

linetype=0 

plot,xcoord2,ycoordl 

plot,xcoord2,ycoordl 

linetype=3 

plot,xcoord2,ycoord2 

linetype=2 

plot,xcoord2,ycoord3 

linetype=5 

plot,xcoord2,ycoord4 

linetype=0 

youts , xmin+110 , .995*ycoordl(0) , 

y outs, xmi n+ 110, .995*ycoord2(0)  , 


'Data' 

'NO' 
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(  xmin 


(  xmin 


( xmin 


else : 


xyout 
xyout 
end 

eq  2200 
!  line 
oplot 
oplot 
!  line 
oplot 
!  line 
oplot 
1  line 
oplot 
!  line 
oplot 
1  line 
xyout 
xyout 
xyout 
xyout 
xyout 
end 

eq  2600 
!  line 
oplot 
oplot 
1  line 
oplot 
!  line 
oplot 
1  line 
oplot 
!  line 
oplot 
!  line 
xyout 
xyout 
xyout 
xyout 
xyout 
end 

eq  3000 
I  line 
oplot 
oplot 
!  line 
oplot 
I  line 
oplot 
!  line 
xyout 
xyout 
xyout 
end 

dummy=0 


s, xmin+110, . 9  9  5*ycoord3 ( 0) , 'LBH' 

s,xmin+110, . 995*ycoord4 ( 0 ) , 'Nil  doublet  2143A' 


)  and  (xmax  eq  2600):  begin 

type=0 

,xcoord2,ycoordl 

,xcoord2, ycoordl 

type=3 

,xcoord2,ycoord2 

type=2 

,xcoord2,ycoord3 

type=l 

,xcoord2,ycoord4 

type=5 

,xcoord2,ycoord5 

type=0 

s , xmin+110 , 

s , xmin+110 , 

s , xmin+110 , 


995*ycoordl( 0) , 
995*ycoord2( 0 ) , 
995*ycoord3( 0) , 
s,xmin+110,.995*ycoord4(0), 
s, xmin+110, .995*ycoord5(0), 


'Data' 

'NO' 

'LBH' 

'VK' 

'Oil  2470A' 


)  and 
type  = 
,  xcoo 
,  xcoo 
type  = 
,  xcoo 
type  = 
,  xcoo 
type  = 
,  xcoo 
type  = 
,  xcoo 
type  = 
s ,  xmi 
s ,  xmi 
s ,  xmi 
s ,  xmi 
s ,  xmi 


(xmax  eq  3000):  begin 


0 

rd 

rd 

3 

rd 

1 

rd 

5 

rd 

4 

rd 

0 

n+ 

n+ 

n+ 

n+ 

n+ 


3 , ycoordl 
3 , ycoordl 

3 ,ycoord2 

3 ,ycoord3 

3 ,ycoord4 

3 , ycoord5 

210, .995*ycoordl(0) 
210, .995*ycoord2(0) 
210, .995*ycoord3(0) 
210, .995*ycoord4(0) 
210, .995*ycoord5(0) 


'Data' 

'NO' 

'VK' 

'01  2972A' 

'OIII  2853A' 


)  and  (xmax  eq  3400):  begin 
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,xcoord2, ycoordl 
,xcoord2, ycoordl 
type=3 

,xcoord2,ycoord2 
type=l 

,xcoord2,ycoord3 
type=0 

s, xmin+110, . 99  5*ycoordl ( 0) , 'Data' 
s. xmin+110, . 99  5*ycoord2 ( 0 ) , 'NO' 
xmin+110, . 99 5 *y coord 3 ( 0 ) , ' VK' 


endcase 


;  c  1  o  s  e  ,  1  3 
close , 12 
close , 11 
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close , 10 
close ,  9 
close ,  8 
close ,  7 
close ,  6 
close ,  5 
close  ,  4 
close ,  3 
close ,  2 
close ,  1 

return 
end 
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